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ABSTRACT 

We discuss the methods employed to photometrically calibrate the data acquired by the Low Frequency Instrument on Planck. Our calibration 
is based on the Solar Dipole, caused by motion of the Solar System with respect to the CMB rest frame, which provides a signal of a few mK 
with the same spectrum as the CMB anisotrophies and is visible throughout the mission. In this data release we rely on the characterization of the 
Solar Dipole as measured by WMAP. We also present preliminary results on the study of the Orbital Dipole, caused by the motion of the Planck 
spacecraft, which agree with the WMAP value of the Solar System speed to 0.2 %. 

We compute the calibration constant for each radiometer roughly once per hour, in order to keep track of changes in the detectors' gain. Since 
non-idealities in the optical response of the beams proved to be important, we implemented a fast convolution algorithm which considers the full 
beam response in estimating the signal generated by the dipole. Moreover, in order to further reduce the impact of residual systematics due to 
sidelobes, we estimated time variations in the calibration constant of the 30 GHz radiometers (the ones with the most important sidelobes) using 
the signal of an internal reference load at 4 K instead of the CMB dipole. 

We estimated the accuracy of the LFI calibration following two strategies: (1) we have run a set of simulations to assess the impact of statistical 
errors and systematic effects in the instrument and in the calibration procedure, and (2) we have performed a number of internal consistency checks 
on the data. Errors in the calibration of this Planck/LFl data release are expected to be about 0.6 % at 44 and 70 GHz, and 0.8 % at 30 GHz. 

Key words, cosmic microwave background - instrumentation: polarimeters - methods: data analysis 



1. Introduction 

This paper, one of a set associated with the 2013 release of 
data from the Planck 1 mission (Planck Collaboration I 2013), 



* Corresponding author: Maurizio Tomasi tomasi@lambrate . 
inaf . it 

1 Planck (http://www.esa.int/Planck) is a project of the 
European Space Agency (ESA) with instruments provided by two sci- 
entific consortia funded by ESA member states (in particular the lead 



describes the techniques we employed to calibrate the volt- 
ages measured by the LFI radiometers into a set of thermody- 
namic temperatures (photometric calibration). We also discuss 
the quality of our calibration in terms of the required accuracy 
needed to achieve Planck's final science goals. This paper is part 
of a larger set of articles (Planck Collaboration II 2013; Planck 



countries France and Italy), with contributions from NASA (USA) and 
telescope reflectors provided by a collaboration between ESA and a sci- 
entific consortium led and funded by Denmark. 
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Collaboration III 2013; Planck Collaboration IV 2013) which 
explain the methodology used to produce maps from raw LFI 
data which have been issued in the Planck 2013 data release. 

By "calibration" here we mean the process that converts 
each digital sample produced by one of the Analogue-to-Digital 
Converters (ADC) into a thermodynamic temperature. The con- 
version can be modelled by the following equation, assuming 
perfect linearity: 



v - G(r s k y + r n0 j se ). 



(i) 



This relates the voltage V measured by the ADC with the temper- 
ature r s k y , obtained by convolving the sky temperature with the 
beam response of the instrument at a given time, and the noise 
temperature r noise of the radiometer. We are primarily interested 
in K = G \ as the purpose of the calibration is to convert V 
back into a temperature. Since a number of environmental fac- 
tors influence the value of K, we expect it to change with time. 
The output of a calibration procedure is therefore a time series of 
values Ki, which sample the unknown function K(t) at a reason- 
able frequency (i.e., higher than the frequency of the expected 
fluctuations in K) and which allow us to reconstruct the value 
of r s k y in Eq. 1 with good accuracy (for LFI, this accuracy is 
between 0.1 and 0.5%). 

Earlier accounts of the calibration procedure for LFI were 
given by Villa et al. (2010) and Mennella et al. (2010), which 
present the results of the LFI on-ground calibration campaign, 
and by Mennella et al. (2011), which describes the LFI cali- 
bration procedure used for producing the Planck Early Results 
(Planck Collaboration I 201 1). 

The structure of this paper is the following. In Sect. 2 we 
introduce a number of important ideas that are going to be 
used in this work, namely the time scales of variations in K 
(Sect. 2.1), the treatment of beam sidelobes in the calibration 
and their impact on subsequent analyses of LFI's calibrated data 
(Sect. 2.2), and an updated list of colour corrections (Sect 2.3). 
Then, in Sect. 3, we explain the methods we have developed to 
calibrate the data acquired by the LFI radiometers. We discuss 
the type of systematic effects affecting the calibration procedure 
in Sect. 4, and we estimate the accuracy of our calibration in 
Sect. 5. Finally, in Sect. 6 we summarize our results and propose 
a number of improvements to be implemented for the future re- 
leases of LFI data. 

2. Calibration philosophy and conventions 

2.1. Time scale of variations in the gain 

In this section we establish the time scale over which we expect 
significant variations in the gain of the LFI radiometers. This 
quantity drives the design of the calibration algorithms we then 
discuss in the next sections. 

Changes in the gain of the LFI radiometers are mainly trig- 
gered by changes in the thermal environment of the LFI in- 
strument, particularly in the front-end and back-end modules 
(Bersanelli et al. 2010). The time scale of gain changes can 
therefore be estimated either by considering the rate of change 
of temperature sampled near the radiometer amplifiers, or by us- 
ing the radiometer to continuously measure the temperature of 
a load kept at a stable temperature. The latter solution is viable 
for LFI, because each radiometer continuously observes a sta- 
ble 4 K load mounted on the external shield of HFI (Valenziano 
et al. 2009; Lamarre et al. 2010); the temperature of each 4K 
load varies by only a few pK per year, and the 4 K signal enter- 
ing LFI radiometers goes through the same chain as the signal 



coming from the telescope. It can therefore be used to assess the 
rate of change in the gain with good accuracy. 

There are a few phenomena that can affect variations in the 
temperature near the amplifiers: 

1 . Fluctuations in the temperature of the 20 K sorption cooler 
cold end (Planck Collaboration II 201 1). The cycle period of 
each cooler's bed is of the order of tens of minutes, but the 
induced change of temperature in the focal plane is minimal 
(less than 0.01 %). We must move to much longer time scales 
(i.e., weeks) to see significant variations in the temperature. 

2. Fluctuations in the temperature of the warm back-end of the 
radiometers. The biggest variations we experienced during 
the mission are due to the continuous turning on/off of the 
transponder antenna (used to send data to Earth) early in the 
mission, which followed a duty cycle of 24 hours and in- 
duced 1 % fluctuations in the temperature near the back-end 
amplifiers. 

We do not expect significant variations in the gain of the ra- 
diometers on timescales shorter than these, i.e., tens of hours. To 
justify this, we can estimate the rate of change of a temperature 
or output voltage /(f) by means of the following parameter 2 : 



fit) 



(2) 



which has the unit of time and quantifies the typical time re- 
quired to induce a change of level s in / at time t. Using s = 0.01 
(i.e., we are looking for 1 % changes), we find the following 
timescales: 

1 . Fluctuations in the temperature Tf p of the focal plane happen 
on timescales of the order of weeks. 

2. The time scale for fluctuations in the temperature 7b e of the 
back-end are faster during the first survey, as (T Tbe ) is of the 
order of tens of hours. After the first survey 3 the transponder 
was left on continuously, and this value increases to roughly 
one week. 

3. Results similar to those for the back-end are found when 
considering the total-power voltage V re f , which measures the 
temperature of the 4K loads (i.e., setting / equal to V re f in- 
stead of Tfp or 7be). 

These results motivate the need to re-calibrate each radiometer 
more than once per day. The most natural length of time for 
Planck is the duration of one pointing period, i.e., the interval 
between two consecutive repointings of the spacecraft, which 
happens roughly once per hour (Dupac & Tauber 2005) and is 
short enough to detect any significant change in the gains. At 
the same time, the interval is long enough to sample the dipole 
signal (our main calibration source) with good signal-to-noise 
ratio, since during one pointing period the telescope scans the 
same circle in the sky tens of times. 

2.2. Handling beam efficiency 

As described by Bersanelli et al. (1997), Cappellini et al. (2003), 
and Mennella et al. (2011), the calibration of LFI is referenced 



2 If the quantity f{t) is affected by statistical noise at high frequen- 
cies, as it is the case for all the quantities considered here, it is necessary 
to apply some kind of low pass filter to it before applying Eq. 2 in order 
to obtain meaningful results. 

3 The number of sky surveys in this Planck data release is 2.5. 
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to the dipole signal, which is a nearly ideal calibrating source. 
For this release, our reference dipole D re f is the combination of 
the solar dipole as given by the WMAP values (Jarosik et al. 
2011), and the orbital dipole, derived from the known velocity 
of the Planck spacecraft relative to the barycentre of the Solar 
System. We now describe the calibration procedure for LFI and, 
in particular, we discuss how non-ideal beams affect the process. 

The output signal (in arbitrary units) measured by a given 
Planck detector pointing at direction x can be written as: 



Voutto = G(B * (r sky + D))(x) + M, 



(3) 



where B is the angular response (beam), normalized so that 
J 4 ^BdD. = 1, D is the dipole (including the solar and orbital 
terms, as well as their associated kinematic quadrupoles), M is 
the monopole (including instrumental offsets), and r s k y is the 
remaining sky brightness temperature (CMB and foregrounds). 
We use the following notation for convolution: 



(A * B)(8, 4>) 



f A(0',cf>' 



)B(8-&,<t>-<p')dQ!. (4) 



Calibration is performed by fitting the measured dipole- 
modulated signal V out (or, in principle, the difference between 
two points in a scan circle) to a model of the beam-convolved 
dipole B mode i * D ref : 



G - G- 



B*D 



B 



model 



* D 



(5) 



ref 



A common approach is to approximate the beam with a pencil 
beam, B mode i ~ Bpentii> as adopted by HFI (Planck Collaboration 
VIII 2013) and WMAP (Jarosik et al. 2007). For LFI we have 
computed the full An beams for all detectors and we developed a 
fast convolution routine to estimate B mod ei * A-ef • This was moti- 
vated by the wish to fully control dipole-coupling to the side- 
lobes, particularly at 30 GHz where sidelobes are larger. The 
convolution of a generic beam with a dipole produces a smear- 
ing effect (due to the fact that not all the power is in the main 
beam) and a slight tilt in the dipole axis (resulting from asym- 
metries in the beam). As shown in detail in Appendix A, these 
effects can be quantified by the length and direction of vector S, 
respectively. 

We verified a posteriori that the net effect of the convolu- 
tion does not produce a significant improvement over a pencil 
beam model at the present stage of the analysis (besides the tilt 
effect, which is < 0.01 %). For this data release, therefore, in the 
convolution routine we rescaled the length of vector S to unity, 
which is effectively equivalent to modelling the beam as a pencil 
beam. The advantage of this approach is that it matches the con- 
vention on the normalization of the beam that has been assumed 
in the calibration of HFI data. In the following description of 
map calibration and associated uncertainties, therefore, we will 
assume B mode i ~ ^pencil- Full use of the 4n convolver, including 
integrated frequency-dependent sidelobes within the LFI radio- 
metric bands, will be applied for later analyses, including polar- 
ization calibration. 

We can write the true beam B as the sum of two terms, B = 
fimain + ^side, where B ma ; n represents the contribution of what 
we define as the main beam (defined as the portion of the beam 
within 5° from the beam centre), and B s i de represents the much 
smaller remaining part (the sidelobes). Equation 5 becomes: 



^ „£main * D + B side * D 

Li ~ (j . 

^pencil * D 



(6) 



The convolution of the main beam with the signal from the sky, 
T - D + r s k y , is nearly identical to an ideal pencil beam except 
that only a fraction 4 1 - / side of the antenna gain is contained in 
the main beam: 

fimain * T « (1 - / si de)Bpencil * T « (1 - f sidt )T. (7) 

So we have: 

,(1 - /side)#pencil * D + B side * D 

(8) 



G - G- 



where 



^pencil * D 
= G(l-/ side Xl+fe), 

_ side * D 

<P °~ B main *D' 



(9) 



Consider now a properly calibrated timeline where the dipole 
and the monopole terms in Eq. 3 have been removed. We can 
write the relationship between the measured sky map T s k y and 
the true sky map r s k y as: 

_ Voul _ £main * ?sky + #side * T sky 
Sky " G = ' d-/side)(l+^) 

Solving for the true temperature map r s k y we find: 



T s ky ~ f s k y (l - 0sky + <Pd), 



where 



0sky = 



sky / sky 



^main * -^sky 



' sky 



(ID 



(12) 



A precise evaluation of the correction terms in Eq. 1 1 would 
require detailed simulations, which have not yet been performed. 
An approximate estimate of <p D can be obtained by taking av- 
erage values of the convolutions involved in its definition. We 
have computed the convolution of the sidelobes with the dipole 
from sample GRASP full beams, leading to <p D as 0.15%. The 
term s k y requires a full convolution of a sky model with the full 
beam, and it is frequency dependent. We performed a simulation 
for a 70 GHz channel and found values for (f> sky ranging between 
0.05 % to 0.2 % throughout a full survey. In conclusion, the cor- 
rection terms s k y and <po are of the same order, within 0.2 %, 
and they tend to cancel out in Eq. 11. Given that the relative un- 
certainties on both terms are large, we do not correct for them in 
the data and estimate a residual uncertainty of 0.2 % in the gain. 
This uncertainty is included in our overall estimated calibration 
uncertainty. 

2.3. Colour corrections 

The raw differential signal measured by a Planck radiometer can 
be written as: 



V = G 1 . >•)'/,,,( ridr = ^- j 'g(v)I(v)A 2 dv 



1 



(13) 



T(v)/(v)dv, 



where G is the overall gain, g(v) is the bandpass, r RJ (v) is the dif- 
ferential Rayleigh-Jeans brightness temperature averaged over 



4 The quantity 1 - / S M e is approximately equal to ||S||. However, as 
explained in Appendix A, the vector S models the coupling of the full 
beam (including the main beam, which causes a tiny smearing effect 
because of its finite FWHM) with the dipole and only uses the first 
directional moments of the beam shape. 
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Table 1. Multiplicative colour corrections C(a) for individual LFI Radiometer Chain Assemblies and for the band average maps. 



RCA Spectral index a 
-2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 

LFI-18 0.948 0.961 0.972 0.981 0.988 0.994 0.997 0.998 0.997 0.995 0.990 0.983 0.975 

LFI-19 0.856 0.878 0.899 0.919 0.939 0.957 0.975 0.991 1.006 1.020 1.032 1.043 1.053 

LFI-20 0.889 0.908 0.925 0.941 0.956 0.970 0.983 0.994 1.003 1.011 1.018 1.023 1.027 

LFI-21 0.917 0.933 0.947 0.960 0.971 0.981 0.989 0.996 1.001 1.004 1.006 1.006 1.004 

LFI-22 1.024 1.026 1.027 1.026 1.023 1.018 1.011 1.003 0.993 0.982 0.969 0.955 0.940 

LFI-23 0.985 0.991 0.996 0.999 1.001 1.002 1.002 1.000 0.997 0.993 0.988 0.982 0.975 

70 GHz 0.938 0.951 0.963 0.973 0.982 0.988 0.994 0.997 0.999 0.999 0.998 0.995 0.991 

LFI-24 0.978 0.984 0.988 0.993 0.996 0.998 0.999 1.000 0.999 0.998 0.996 0.993 0.989 

LFI-25 0.967 0.974 0.980 0.985 0.990 0.994 0.996 0.999 1.000 1.000 1.000 0.999 0.997 

LFI-26 0.957 0.966 0.973 0.980 0.985 0.990 0.995 0.998 1.000 1.001 1.002 1.002 1.000 

44 GHz 0.968 0.975 0.981 0.986 0.990 0.994 0.997 0.999 1.000 1.000 0.999 0.998 0.995 

LFI-27 0.948 0.959 0.969 0.978 0.985 0.991 0.995 0.998 1.000 1.000 0.998 0.995 0.991 

LFI-28 0.946 0.958 0.968 0.977 0.985 0.991 0.996 0.998 1.000 0.999 0.997 0.993 0.988 

30 GHz 0.947 0.959 0.969 0.977 0.985 0.991 0.995 0.998 1.000 1.000 0.998 0.994 0.989 



the beam, and G' is defined such that 2k B G'r(v) = Gg(v)A 2 . The 
constant G' would be the gain if bandpasses were defined via 
the transmission coefficient t(v) oc g(v)A 2 instead of g(v), as HFI 
does (Planck Collaboration VIII 2013). 

Since we calibrate using the CMB dipole, the calibration sig- 
nal measured in any given pointing period is: 



d r 

K=°} 



g(v)D T] AT (v)dv, 



(14) 



where D is the dipole amplitude in thermodynamic units and 
i]at(v) is the conversion factor from CMB temperature to RJ 
temperature. Hence the calibrated sky map temperature is: 



T = KV 



/ g(v)T R3 (v)dv 



(15) 



If we are observing pure CMB fluctuations, then T RS (v) = 
Ar r/ AT (v), and hence we have T = AT as expected. If we are 
observing foreground emission with intensity power-law index 
a, then T RJ (v) = T (v/v r~ 2 , and 

Jg(v)(v/vor' 2 dv 

/ g(y)j]AT(y)dy 

The LFI colour correction is defined as: 

/ g(v)77Ar(v) dv 



T = T< 



(16) 



C(a) 



?7Ar(vo) / g(y)(y/vo) a 2 dv 
Hence the colour-corrected temperature is: 
C(a)T = T /t] AT (vo). 



(17) 



(18) 



This gives the thermodynamic brightness differential tempera- 
ture (units "Kcmb") at the reference frequency vo. To get to the 
standard (Rayleigh-Jeans) brightness temperature at the same 
frequency we use rj AT {yo) (see Eq. 14): 



r [K RJ ] = f[K C MB]T]AT(y )C(a). 



(19) 



For the present series of Planck papers, the reference frequencies 
v are defined as precisely 28.4, 44.1, and 70.4 GHz for the three 
LFI bands; in the Planck Early Release papers, slightly different 
values were used. 



Values for C(a) for the expected range of foreground spec- 
tral indices are listed in Table 1 . Note that this definition of C is 
inverted relative to the colour correction C(a) used in the Planck 
Early Release. We note additionally that C(a) quoted here are 
not exactly equal to 1/C(ar) from Mennella et al. (2011) and 
Zacchei et al. (2011), due both to the change in reference fre- 
quency and also to a small error in the earlier estimates. 

Colour corrections at intermediate spectral indices may be 
derived accurately from a quadratic fit to the values in Table 1 . 
In addition, the data release includes the Unit conversion and 
Colour Correction (UcCC) IDL package (Planck Collaboration 
IX 2013) which calculates colour corrections and unit conver- 
sions using the band-averaged bandpass information stored in 
the Reduced Instrument MOdel (RIMO) file also included in the 
data release. 

From Eq. 17 it appears that the absolute scaling of the band- 
pass is irrelevant for the colour correction, although Eq. 13-15 
require that: 



g(y) r]Ar(y) dv = 1, 



(20) 



for consistency with the definition of gain, G, and calibration 
factor, K, elsewhere in this paper. In practice, all the colour cor- 
rections listed in Table 1 are derived from averages across two 
or more bandpasses: values for individual Radiometric Chain 
Assemblies (RCA) require averaging the main and side radiome- 
ters in each RCA, and the response of each radiometer is the 
average of the two independent detectors (Zonca et al. 2009). 
For consistency, the weighting for these averages must duplicate 
the procedure used to average the data, as described in Planck 
Collaboration II (2013). We recall here the procedure: (i) cali- 
brate the individual data streams to give equal response to the 
CMB; (ii) combine the raw data for the detectors for each ra- 
diometer using fixed weights based on the inverse variance mea- 
sured early in the mission, incorporating fixed calibration val- 
ues; (iii) recalibrate the data using time-dependent factors as de- 
scribed in the present paper; (iv) combine the two radiometers in 
each RCA with equal weights to minimize polarization leakage; 
and (v) combine the data from individual RCAs with inverse- 
variance weights (fixed for the whole mission). When combin- 
ing the bandpasses, the initial calibration step is equivalent to 
normalizing each component bandpass according to Eq. 20. The 
band-averaged bandpass stored in the RIMO, however, normal- 
ized the bandpasses using J g(v) dv = 1. The difference is minor, 
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since ij&t varies by only a very small amount within any of the 
LFI bands, but it accounts for a small (< 0. 1 %) difference be- 
tween the C(a) values listed here and those derivable from the 
RIMO bandpasses. 

The main source of uncertainty in the values of C(a) are the 
uncertainties in the bandpass (Zonca et al. 2009), which were not 
well quantified; flight data now allow an indirect estimate, as fol- 
lows. The two radiometers in each RCA, known as the main- and 
side-arms, are sensitive to orthogonal polarizations. The band- 
passes for the two arms differ, leading to different colour cor- 
rections. The polarization signal is derived from the difference 
of the calibrated signals from the two arms (Leahy et al. 2010); 
unpolarized foreground emission does not precisely cancel due 
to the diffential colour factors. This is the "bandpass leakage" 
effect. This leakage can be estimated from the flight data, as de- 
scribed by Planck Collaboration II (2013), but it can also be es- 
timated from the pre-launch bandpass models that we use here 
to calculate C(a). To a good approximation, we can write the 
leakage factor as 

(Cs(a)-Cu(a))/2 = (J3-J3 CM B)a, (21) 

where the a-factors depend solely on the bandpass profile, B = 
a - 2 is the temperature spectral index, and /?cmb is the in-band 
spectral index of the CMB, which we can take as zero for present 
purposes. The flight measurements of the a-factors are demon- 
strably more accurate than the pre-launch ("QUCS") estimates, 
and so 

08 - y8cMB)(a Q ucs - flfl lg ht) ~ (5C s (a) - 6C M (a))/2, (22) 
where 6C is the colour correction error. Statistically 

((5C s (a) - 6C M (a)) 2 ) = ((6C s (a) + SC M (a)) 2 ) = a*, (23) 

where cr c is the colour correction error for one radiometer pair, 
i.e. one RCA. Hence for individual RCAs, <x c ss \B\ <r a . There are 
too few pairs for this approach to give an accurate value for cr c , 
especially as the error sources for the bandpasses in each band 
are different; therefore results from the three bands cannot be 
combined. Our rough estimates for cr a are 0.14%, 0.46%, and 
0.5 1 % in the 30, 44, and 70 GHz bands, respectively. For the 
band-averaged maps, the errors in the colour corrections are re- 
duced by V-Wrca, giving overall rough uncertainties of 0. 1 1>8|, %, 
O.3|j0| %, and 0.2\/3\ %, respectively. Given typical values of B of 
around -2 to -3, this gives uncertainties in the colour correc- 
tions of a few tenths of a percent at 30 GHz, and 0.5-1 % at 44 
and 70 GHz. 

Progress towards in-flight calibration of the colour correc- 
tions is discussed further in Sect. 5.4.1. 

3. Calibration techniques for LFI 

There is not a unique way to produce an estimate of the calibra- 
tion constant K (see Eq. 1). In picking the method to use, one 
has to consider a number of elements: 

1. the calibration should be as accurate as possible; 

2. calibrated data should also be self-consistent; for example, 
with the same beam and the same pointing direction and ori- 
entation in the sky, the value of T^ y at any frequency should 
not depend on when the measurement was done, nor on the 
detector (under the hypothesis that everything else, e.g., the 
bandpass, can be assumed to be the same); 

3. the estimated shape of K(t) should be motivated by a physi- 
cally meaningful model of the radiometer. 



To help in the classification of the calibration methods discussed 
in this work, we write: 

K(t) = K (l+t(t)% (24) 

thus decomposing K into a constant term Kq and a unitless time- 
varying quantity £(f). Typically £(f) varies by a few percent per 
year. This allows us to categorize calibration methods into the 
following families: (1) absolute methods produce an estimate 
for Kq\ (2) relative methods estimate how £ changes with time; 
and (3) some methods are able to estimate K(f) directly and can 
therefore be considered both absolute and relative. 

We implemented two calibration schemes for the LFI ra- 
diometers: 

- The OSG dipole calibration is an improved version of the 
pipeline used to calibrate the data used in the Planck Early 
and Intermediate data releases. It is based on the OSG algo- 
rithm ("Optimal Search for Gains") and relies on the signal 
of the solar and orbital dipoles as observed by the spacecraft. 
This method is only weakly affected by non-idealities in the 
radiometers (e.g., non-linearities in the ADCs, see Sect. 4.3), 
but optical effects (e.g., sidelobes, see Sect. 4.1) can induce 
systematic errors in the reconstruction of the gain. 

- In the 4 K calibration we use the dipole only to fix the ab- 
solute level of calibration (Kq in Eq. 24), but we estimate 
gain changes (the £ term) using the 4 K total-power output 
of the radiometers. Unlike the OSG dipole calibration, this 
scheme is sensitive to ADC non-linearities, but it provides 
an estimate for the gain K which is independent 5 of optical 
effects. 

The chief reason why we decided to use two different calibra- 
tion schemes was to improve the self-consistency of the maps. 
We were not able to derive a calibration method that was robust 
against both radiometric non-idealities and optical effects in time 
for this data release. Since optical effects are most significant in 
the 30 GHz channels, we decided to apply the 4 K calibration to 
these radiometers instead of the other one (which was the base- 
line for the Early and Intermediate PlanckfLFl papers). We pro- 
vide a summary of the calibration methods used for each LFI 
frequency in Table 2. In Sect. 4.1 we will discuss the level of 
consistency between the two calibration methods in the context 
of the treatment of optical systematics. 

The outline of this section is as follows. The OSG dipole 
calibration and the 4 K calibration are explained in Sect. 3.1 and 
3.2. In Sect. 3.3 we present the current status of our efforts to 
use the orbital dipole signal for calibrating data. Finally, Sect. 3.4 
deals with the task of setting the zero-level of Planck's LFI maps. 

3.1. OSG dipole calibration 

The combination of the motions of the spacecraft, the Earth, and 
the Solar System with respect to the CMB produce a dipole sig- 
nal with an amplitude of (3.355 + 0.008) mK (Hinshaw et al. 
2009) on the full sky. Like COBE (Kogut et al. 1996) and WMAP 
(Jarosik et al. 2011), we use this signal to estimate the value of 
K in Eq. 1. We indicate the gain estimates produced using this 
method with K& v . 

The basic steps of this calibration procedure are as follows: 

5 Note however that, even if the calibration constants A", are not af- 
fected by uncertainties in beam shapes, this is not necessarily true for 
the calibrated data. After having applied the calibration factors, we must 
remove the dipole signal convolved with the beam from the data, and 
therefore any uncertainty in the beam shape will lead to a systematic 
error in the calibrated data, even if the 4 K calibration is used. 
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Fig. 1. Variation in time of a few quantities relevant for calibration, for radiometer LFI21M (70 GHz, left) and LFI27M (30 GHz, 
right). Vertical dashed lines mark boundaries between complete sky surveys. All temperatures are thermodynamic. Panel A: 
Calibration constant K estimated using the expected amplitude of the CMB dipole. Note that the uncertainty associated with the 
estimate changes with time, according to the amplitude of the dipole as seen in each ring. Panel B: Expected amplitude of the dipole. 
The shape of the curve depends on the scanning strategy of Planck (see also Fig. 3), and it is strongly correlated with the uncertainty 
in the gain constant (see panel A). The two dots indicate at which time the data used to produce the plots in Fig. 4 were taken. Panel 
C: 4 K total-power output voltages V. Note that the noise level for V is much more stable than it is for K , as it does not depend 
on the dipole amplitude. (All the large jumps seen in this figure occurred because of some change in the operational state of the 
spacecraft.) Panel D: The calibration constants K used to actually calibrate the data for this Planck data release are either derived 
by applying a smoothing filter to the raw gains in panel A (70 and 44 GHz, e.g., LFI21M, left), as described in Sect. 3.1, or through 
Eq. 30 and the 4 K total-power voltages in panel C (30 GHz, e.g., LFI27M, right) 

Table 2. Methods used to calibrate the radiometers. 



Frequency 


Absolute cal. 


Relative cal. 


Primary systematics 


Relevant paper sections 


30 GHz 


CMB dipole 


4K 


ADC non-linearities, variations 


3.1, 3.2, 4.3 








of r noise with time 




44 GHz 


CMB dipole 


CMB dipole 


Optical systematics, bandpass 


3.1,4.1 








response of the radiometers 




70 GHz 


CMB dipole 


CMB dipole 


Same as for 44 GHz 


3.1,4.1 



1. We combine the velocity of the spacecraft with respect to 
the Sun Vpianck and the velocity of the Sun with respect to 
the CMB Vsun to estimate the amplitude and alignment of the 
dipole in the sky in a given direction x: 



AT(x, t) = T t 



CMB 



1 



7(0(1 -110(011 cos 6{ij) 



1 m) ■ x, 



(25) 

where Tqmb is the temperature of the CMB monopole, vs un 
is the velocity of the Solar System with respect to the CMB 
rest frame, v P i anck is the spacecraft's velocity with respect to 
the Solar System's barycentre, /3 = (vs un + Vpianck)/c, y = 
(1 - B 2 y l/Z , and is the angle between x and /}. Table 3 
reports the numerical values used in the pipeline. 
2. We produce discrete Time Ordered Data (TOD) of the ex- 
pected overall dipole signal (thermodynamic temperature) at 
time tf. 

ATi = AT(x(ti), t-), (26) 



where i ranges from 1 to N ~ 10 5 (the number of sam- 
ples in a pointing period). As we explained in Sect. 2.2, 
the computation of the expected dipole signal takes into ac- 
count the shape of the beams, following a method described 
in Sect. 4.2. 

3. Using pointing information, we project both V, (the voltage 
V, at time i, as used in Eq. 1) and ATj on a HEALPix map 
(Gorski et al. 2005) with N side = 256. Multiple hits on the 
same pixels are averaged in both cases. The result is a pair 
of maps, y™ dp and Ar™ dp , with k being the pixel index 6 . 

4. We use weighted least squares to estimate K in Eq. 1 from 



the correlation between the signal in volt Vf ky with Ar 



sky, 



ymap 



* dlp A7T P + e , 



(27) 



6 Most of the pixels in the maps contain no data, since during one 
pointing period the beam paints a thin circle in the sky. We assume 
hereafter that the index k runs only through the pixels which have been 
hit at least once. 
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Table 3. Parameters used in the model of the dipole signal. 



Parameter Value Source 

CMB monopole 
Jcmb 2.725 ± 0.002 K Mather et al. (1999) 

Solar speed" 

/ 263!99 ± 0? 14 Jarosik et al. (201 1) 

b 48!26 ± 0?03 Jarosik et al. (201 1) 

v Sun 369.0 ±0.9 km s~> Jarosik et al. (201 1) 

Velocity of the spacecraft 6 

<v P i anck > 30.0 ±0.4 km s" 1 

minv P i anck . . . 29.39 km s -1 
maxvpi anck . . . 30.60 km s _I 



8 Relative to the CMB rest frame. 

b The values reported here are representative of the overall speed of the 
spacecraft. The calibration code uses the full timestream of velocity 
components (y x , v y , V z ), sampled once every minute. 




Fig. 2. Mask used in the application of Eq. 27 to 30 GHz data 
(Ecliptic coordinates). The mask hides point sources and the 
strong emission of the Galactic plane. What is left is that part 
of the sky where the only significant emissions are the CMB 
dipole and the CMB itself. Similar masks have been used for 44 
and 70 GHz data. 



where ^ip and e are the parameters used in the fit. Each sam- 
ple k is weighted according to the number of hits per pixel. 
In computing the fit, we use a frequency-dependent mask to 
avoid those pixels where a strong non-Gaussian signal other 
than the dipole is expected, i.e., point sources and the Galaxy 
(see Fig. 2). 

5. The fit described in the previous step suffers from the pres- 
ence of noise in the sky signal (mainly due to the CMB and to 
Galactic emission that was not properly masked). Therefore 
we apply an iterative algorithm, named Mademoiselle that 
uses destriping techniques to iteratively improve the map 
(Planck Collaboration II 2013) by constraining multiple 
crossings of the same pixel to have the same signal by a 
(small) shift in their gains. 

6. The gains produced by this procedure need to be further pro- 
cessed in order to reduce the statistical noise. We applied an 
adaptive smoothing filter based on wavelets which smooths 
more around dipole minima and less around maxima (see 
Planck Collaboration II 2013). 

7. Once a set of K, gains (one per pointing period) is produced, 
the pipeline calibrates the timelines and subtracts the dipole 
signal, so that in the calibrated timestreams the Galactic sig- 
nal and the CMB fluctuations are the only relevant astrophys- 
ical components. 




50 100 150 200 250 300 350 
North phase angle [deg] 



Fig. 3. Top: dipole signal due to the motion of the Solar System 
only (in Ecliptic coordinates). The points observed by LFI27M 
during pointing periods 4800 (250 days after launch) and 6636 
(300 days) are shown with black and blue lines, respectively. (In 
Fig. 1, panel B on the right, we indicate the times when these two 
periods occurred with black dots.) Bottom: dipole signal along 
the same two pointing periods as a function of the North phase 
angle. Note the difference in the amplitude of the two sinusoidal 
waves. Such differences are due to the Planck's scanning strat- 
egy and determine the statistical error in the estimation of the 
calibration constant K = G l (Eq. 1). Refer also to Fig. 4 for 
details about how K is computed for these two pointing periods. 



Details of the code implementation are given in Planck 
Collaboration II (2013). 

The accuracy in gain reconstruction depends critically on the 
orientation of the spacecraft with respect to the dipole axis, i.e. 
the value of the product (vs un + Vpi an ck) ■ x. We speak of a dipole 
minimum when the spacecraft's orientation is such that the value 
of the scalar product reaches a minimum, and a dipole maximum 
when reaches a maximum. Figure 1 shows this idea for two LFI 
radiometers, LFI21M (left) and LFI27M (right): in panel A we 
show the values of K,a v as a function of time, while in panel B 
we show the expected amplitude of the dipole Ar<ji P as seen in 
the circle in the sky (ring) which is scanned during each pointing 
period. It is clear that the noise in K& v is mainly due to the vari- 
ation of A7di p . The reason why the peak-to-peak difference of 
the measured dipole signal varies is represented in Fig. 3, which 
shows which parts of the sky are observed by LFI27M during 
two pointing periods, one (4800) near a dipole maximum, and 
another (6636) near a dipole minimum. Finally, Fig. 4 shows the 
data used to compute the weighted linear regression presented in 
item 4 above for the same pointing periods. 
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Fig. 4. The iterative calibration procedure fits the differential 
voltages with the expected dipolar signal (K CMB ), both projected 
on a map, and it calculates the gain according to Eq. 27. As Fig. 3 
shows, depending on the position of the circle with respect to the 
axis of the dipole signal, the range spanned by the signal can be 
either large (dipole maximum, top) or small (dipole minimum, 
bottom). Of course, in the latter case the estimation of the gain 
suffers from a larger statistical error. 



3.2. 4 K calibration 

Gain changes in a LFI radiometer can be modeled using the 
emission of the internal 4 K reference load as measured by the 
radiometer itself (Valenziano et al. 2009). 

We indicate with V ref the output voltage which measures the 
temperature of the 4 K reference load. Equation 1 holds for T re f 
as well, of course. Solving for K = G 1 yields a formula that 



can be used to estimate K from the ratio between the system 
temperature and the 4 K reference voltage: 

. . ^ref + ^iso^sky + T x 

K 4K (t) = 

Vref 



' noise _ ^sys 



Vref' 



(28) 



where e; so is the isolation (for its definition see Villa et al. 2010) 



and the term t^o^sky represents the amount of signal coming 
from the sky and leaking into the reference load signal. We ex- 
press the term V~l as: 



ref 



v ref,0 



1 + 



Vref - V 



ref.O 



V 



ref.O 



1 - 



V re f,0 
Vef - V 



Vef 
Vef 
V re f,0 



,0 / 



(29) 



where we applied a first-order Taylor expansion 7 for V reI around 
the value V re f,o- We substitute it into Eq. 28 to derive the follow- 
ing equation: 



K 4K (t) * 



^iso ^sky^ 



V 



ref.O 



2- 



VreHfV 
V, 



ref.O 



(30) 



(Ideally, variations in T le f, in €i s0 T re f, and in r no i se should be neg- 
ligible; we will quantify their variation later.) In principle, Eq. 30 
would be enough to estimate K at any time, if one knew all the 
terms in the right side of the equation with the desired accuracy 
(e.g., 0.1 %). Since this is not our case, to find the unknown fac- 
tor in the equation we use the estimate of K(t) from the dipole 
(see Sect. 3.1) to perform a weighted linear least squares from 
which we estimate Kq: 



K div (t) = K x\2- 



Vef(0 
V re f,0 



(31) 



(Note that here we use the "raw" values for ^di p (f), i-e., before 
applying the smoothing filter described in Sect. 3.1, point 6.) 
The average value 8 of V re f(f) over the whole set of data to be 
calibrated is used for V re f,o, an d tne weights w, in the linear 
squares fit are proportional to the amplitude of the dipole sig- 
nal Ar dip seen by the radiometer during the z'-th pointing, that is, 
Wi oc Ar^. According to our classification based on Eq. 24, this 
calibration scheme is relative (with £(f) = (V re f,o ~ V re f(f))/Vref ; o) 
but not absolute, as we need K d i v (t) in order to estimate K . 

As an example, see Fig. 1, panel D (right), which shows the 
values of K 4K calculated for LFI27M (a 30 GHz radiometer). 
Unlike the case for LFI21M (panel D, left), which required a 
smoothing filter to be applied to the raw gains in panel A, here 
we use the total-power voltages without any filtering, as their 
error does not depend on the amplitude of the observed dipole in 
the sky, but only on the stability of the 4 K reference load, which 
is better by orders of magnitude. 

7 The error from the first-order approximation in Eq. 29 is of the order 
of 0.05%. 

8 In principle any value for K ref might be used. However, the fact 
that we derived Eq. 30 using a Taylor expansion of K(t) around V rc f,o 
suggest that we use a value of the voltage that is the closest to all the 
values of V(t) we expect during the survey, in order to improve the 
accuracy of the first-order approximation. Depending on what we mean 
by "closest", we might either choose the median or the mean of V(f); 
in this context we found that the difference between these quantities is 
probably negligible, as it is always less than 0.5 %. 
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Fig. 5. Variation of the temperature of the back-end with time, 
as measured by sensor LM207332 (placed near the back-end am- 
plifiers of LFI27M and 27S). The two major thermal events that 
happened during the period of the nominal mission are clearly 
visible here: (1) 258 days after launch the transponder antenna 
switched from a on/off duty cycle of 24 h to being always kept 
on, and (2) 456 days after launch we switched to the other sorp- 
tion cooler used to cool the focal plane. 

When we presented Eq. 30, we asked that variations with 
time of the three terms T re f, ei so r sky , and r no ; se be negligible, 
in order not to induce significant systematics in the determi- 
nation of the gain. We will now quantify how much they are 
expected to change during the mission. (Refer to Sect. 2.1 for 
a general discussion of time scales and LFI calibration.) If we 
write Vref(0- V re f,o as <5V re f, and similarly for 6T lef , 6T sky , 6T noise , 
and 6K, from Eq. 30 we derive the following expression: 

6K(t) _ SV lef + ST Tei + 6 iso 6T sky + 6T noise 

Ko Vref.O T tou0 

withTVo = r ref> o+e iso r sky> o + r noisej o ~ 20 K for the LFI 30 GHz 
radiometers. This equation expresses fluctuations in K in terms 
of fluctuations in V re f and the three terms in the numerator of 
the second fraction on the right side, whose sum should be much 
smaller than r to t,o- In the case of LFI's 30 GHz radiometers, the 
following considerations apply: 

1. The variation of T K f during Planck's nominal mission is less 
than 10 mK, therefore making the term 5T re f/T tot fi negligible 
(around 0.05 %). 

2. The value of ei so is generally less than 0.1 (Villa et al. 2010). 
The biggest contribution to c5r sky comes from the solar dipole 
(max(r sky (f) - r sky0 ) ~ 3.5 mK). Such variations have an 
impact of only about 0.01 % on K^k- 

3. For #r n oise/7tot,o the approximation is not so good, because 
of the dependence of r no ; se on the temperature of the focal 
plane (Terenzi et al. 2009). Two events triggered significant 
changes in the satellite thermal environment. Firstly, starting 
258 days after launch, the transponder of the communica- 
tion system was always kept on, instead of being switched 
on just for transmission, in order to increase thermal stabil- 
ity. Secondly, 456 days after launch, the first sorption cooler 
reached its end-of-life and the second sorption cooler was 
switched on (see Fig. 5). These major changes of the ther- 
mal environment are expected to have an impact on r noise of 
0.5-1 %. This cannot be considered negligible and has to be 
corrected for. 



The simplest correction strategy for this last effect would be 
to apply the 4K calibration model separately for each section 
of the data. However, this technique is more subject to the sys- 
tematic effects known to affect the dipole calibration. Therefore 
a more robust approach is to create a map of each section inde- 
pendently, removing both the cosmological and orbital dipoles, 
then fit for a residual dipole while masking the Galactic emis- 
sion, and finally estimate the correction factor £t K as: 

&k(0 = 1 - — . (33) 

A-/ S k y 

where AT& t (i) is the dipole fitted on each section i and Ar sky 
is the cosmological dipole. For example, if the residual dipole 
is 0.1 % of the cosmological dipole, it means we are over- 
calibrating by 0.1 %, so we need to correct the calibration by that 
factor. Table 4 lists the residual dipole values and the corrections 
applied. 

3.3. Orbital dipole calibration 

The CMB dipole induced by the peculiar velocity of the solar 
system relative to the CMB last scattering surface is a superb tool 
for LFI calibration, thanks to its high signal to noise source for 
all the frequency channels of LFI. Unfortunately it is not a fun- 
damental cosmological parameter and has no predicted ampli- 
tude and direction: therefore, it must be calibrated against some- 
thing else which is known absolutely. The relative velocity of 
the Planck satellite with respect to the solar system barycentre 
is known very precisely and adds a further dipole signal with a 
known amplitude and direction to the data. As the satellite is lo- 
cated at L2, its orbital velocity is similar to the Earth's, that is, ap- 
proximately 30 kms -1 , and the amplitude of the induced dipole 
signal is 250 /iK. This is an order of magnitude smaller than the 
cosmological dipole, and it is too weak to use on a pointing-by- 
pointing basis. Therefore, some form of relative calibration must 
be used to tie together data from each pointing period. 

WMAP 5-year and 7-year data (Hinshaw et al. 2009; Jarosik 
et al. 2011) used a parametrized fit, using detector voltages and 
hot and cold stage temperatures, to provide a relative gain solu- 
tion to which the orbital dipole provided the absolute tempera- 
ture to voltage calibration. The psuedo-correlation radiometers 
of LFI, which use a reference load, are more complicated than 
the WMAP switched differential radiometers. As such, the pa- 
rameterization of the radiometer gain has not reached an usable 
state yet; the 4K calibration (Sect. 3.2) is the first implementa- 
tion. 

A scheme to use the orbital dipole must start with the time- 
ordered-data, since the direction of the induced dipole is in the 
direction the satellite is traveling and therefore is time depen- 
dent. Furthermore the TOI before the calibration module must 
be used as it removes both dipoles from the data. As most of the 
pixels on the sky are observed with two observations separated 
by approximately six months, where the velocity is in opposite 
directions, the orbital dipole signal cancels out to first order on 
the full mission sky maps. The time-ordered data are binned into 
rings for each pointing period, based on the rotation phase angle 
of the satellite. In this way all the relevant information is stored 
in a more compact and accessible form, with the advantage that 
the orbital dipole appears as a near sinusoid variation with the 
same phase in each ring. 

A good relative calibration is required to remove the gain 
differences the between rings. For this we use the gains used 
in the final analysis as these are the best available. For the 70 
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Table 4. Dipole residuals and correction factors applied to the 4 K calibration at 30 GHz in response to major thermal events. 



Time range [days] 


Dipole LFI27 


Dipole LFI28 


Correction LFI27 


Correction LFI28 


91-257 


-17.03 yuK 


-8.91 fiK 


0.58 % 


0.30% 


258-454 


-17.62 yuK 


-13.06 ytiK 


0.60% 


0.44 % 


455-563 


20.66 fiK 


-6.71 //K 


-0.70% 


0.23 % 



Table 5. Preliminary estimates of the CMB dipole in the 44 GHz 
channels. 



Horn 


Arm 


Correction" 


v Su n[kms '] 


I 


b 


24 


M 


1.0005 


367.76 


263?96 


48!54 




S 


1.0000 


368.03 


263!96 


48!50 


25 


M 


1.0006 


367.97 


263!93 


48!45 




S 


0.9967 


369.35 


264! 01 


48!28 


26 


M 


1.0025 


367.44 


263! 87 


48!42 




S 


1.0004 


368.33 


263!92 


48!41 



" The correlation factor of expected orbital dipole to measured dipole. 



and 44 GHz channels this may appear to be potentially circular 
as they are already calibrated on the solar dipole, but tests intro- 
ducing a scaling factor into the relative gains shows that the same 
factor is recovered in the absolute gain. In this way any tension 
between the WMAP dipole used in the calibration and the orbital 
dipole would be revealed as a non-unity correction factor. 

Both solar and orbital dipoles can then be fitted in the time 
line, with the direction and amplitude of the solar dipole left 
free and the amplitude of the orbital dipole providing the final 
calibration. At the moment the uncertainties in the ring-to-ring 
relative calibration and far-sidelobes limit the accuracy of this 
calibration to 0.26 %, 0.16 % and 0.64 % at 70, 44 and 30 GHz, 
respectively, based on variations in the results. As this does not 
offer a significant improve over the WMAP dipole, this analysis 
has therefore been postponed to a subsequent release of the LFI 
data when the characterization of far sidelobes and gain varia- 
tions is more mature. Table 5 reports the results of a preliminary 
analysis of correction factors and solar dipole parameters for the 
44 GHz channels (where the effect of far-sidelobes is least). 

3.4. Setting the zero levels in the maps 

The monopole of Planck maps is unconstrained: adding a con- 
stant to the time-ordered data leaves the likelihood of the map 
unchanged. It is therefore conventional to adjust the monopole of 
the maps in order to attain plausible absolute values of the high 
Galactic latitude diffuse Galactic foreground signal, and plausi- 
ble values of the frequency scaling of this Galactic signal from 
channel to channel. 

For LFI we have implemented two different approaches to 
evaluate the zero level. The first considers a circular region of 
sky centred on (/, b) = (344?47, -77?08) with radius of one de- 
gree 9 , on which the monopole matching is performed. We sub- 
tracted the estimated contribution of the CMB (smoothed to the 
angular resolution of each channel) from each LFI frequency 
map prior to matching diffuse foreground amplitudes. The aver- 
age value of the CMB-subtracted pixels within the above mask 
is adjusted to match following values: 35 //K at 30 GHz, 15 /J.K 
at 44 GHz, and 1 8 //K at 70 GHz. 



9 We chose this region since it is nearly free of bright point sources 
and is well representative of the Galactic diffuse emission. 



The other approach we implemented is similar to the one 
adopted by WMAP and described by Bennett et al. (2003, 2012). 
The procedure is the following: (1) we smooth the maps at 
1° angular resolution, (2) we subtract the estimated CMB sig- 
nal, again smoothed at the some common resolution, and (3) 
we fit the observed variation with Galactic latitude assuming 
a plane-parallel model for the Galactic emission. We consider 
only the behaviour on the southern hemisphere (in the range 
-90° < b < -15°) in order to avoid possible contamination 
from the high-galactic latitude structures present in the Galactic 
emission such as the North Polar Spur. We divide the map into 
stripes of constant latitude 2.5° wide and we evaluate the mean 
level of CMB subtracted signal. Then, this is fitted with a cose- 
cant model of the form T — A esc b + B. The offset is the value 
by which the intercept B is zero. 

The results yielded by the two methods are in agreement. We 
set the zero levels of the maps released by the Planck collabo- 
ration to the estimates obtained using the second method. The 
numerical values are reported in Table 8. 



4. Systematic effects that affect calibration 

Raw gains from the nominal pipeline are sensitive to various ef- 
fects related both to previous data processing steps and to instru- 
ment characteristic. A full treatment of systematics in the cur- 
rent release of LFI data is provided by Planck Collaboration III 
(2013); here we concentrate only on those systematics that have 
a significant impact on the calibration of the radiometers. 

4.1. Impact of the optics on calibration estimation 

The optical response of each radiometer, including the effect 
caused by the presence of the telescope and the baffles, impacts 
the estimation of the calibration constants. Apart from the ef- 
fects due to the coupling of the CMB dipole with the full beam 
response described in Sect. 2.2, an important effect (particu- 
larly significant at 30 GHz) is the fact that the brightness of the 
Galactic plane produces a non-trivial spurious signal which bi- 
ases the fit with the dipole. 

Consequently, for the current release, we have decided to 
calibrate the 30 GHz channels using the 4 K calibration, as de- 
scribed in Sect. 3.2. This procedure avoids the spurious signals 
induced by the sidelobes 10 (see Fig. 6). In a future release, we 
plan to use the fully modelled sidelobe patterns to correct the 
timelines in the iterative calibration sequence, thereby using the 
same calibration scheme for all three frequency bands. We will 
also treat polarization calibration. 

The optical modelling of LFI, including extensive discussion 
of the sidelobes, is covered in Planck Collaboration IV (2013). 



10 Note however that the two methods produce maps that show a re- 
markable agreement: their peak-to-peak difference, if we neglect a few 
outlying pixels, is of the order of a few tens of /iK, thus of the same or- 
der of magnitude as the r.m.s. of the survey difference maps discussed 
in Sect. 5.3.2. 
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Table 6. Dipole deflection and correction due to sidelobes. 



Radiometer 
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Fig. 6. Difference between the nominal 30 GHz LFI map (cali- 
brated using the 4K method, see Sect. 3.2) and the same map 
produced using the OSG calibration (Sect. 3.1). The residual 
dipole visible in the map is clearly a sign that one of the two 
calibrations (OSG) is not able to properly fit the dipole. A few 
stripes are visible as well: they are caused by Galactic straylight 
affecting the fit with the dipole, a problem which mainly affects 
the OSG calibration. Similar maps have been used to calibrate 
44 GHz and 70 GHz data. 



4.2. Efficiently convolving the dipole with a realistic beam 

In Sect. 2.2 we mentioned that the procedure is based on a convo- 
lution of a reference dipole with the beam shape. This procedure 
can be time-consuming, as it must be done once for each sample 
measured by each radiometer. Fortunately, the spatial symmetry 
of the dipole signal allows us to derive a fast and elegant method 
to properly take into account the shape of the full beams when 
using the dipole to calibrate the data. We describe the mathe- 
matical details of the model in Appendix A; here we provide a 
physical interpretation. 

The model applies a convolution with the full beam to 
Eq. 25, and it takes advantage of the high level of symmetry 
in the dipole to greatly simplify the calculation. The result of the 
convolution is still a dipole signal, but its direction and amplitude 
are slightly different from the dipole before the convolution, for 
the following reasons: 

1 . the amplitude changes because of all the power entering the 
beam from directions other than the beam axis; 

2. the direction changes because of asymmetries in the beam 
shape. 

Such changes are characterized by a constant vector 5, which 
only depends on the beam shape. The maximum difference be- 
tween the amplitudes of the convolved dipole and of the dipole 
in the sky is: 



Ar sc - 7cmb 



VSu 



(i - \s\). 



(34) 



The direction of S in the beam reference frame represents the 
deviation of the effective axis of the dipole signal with respect 
to the axis of the dipole in the sky. A convenient measure of the 
overall deflection is represented by the deflection angle 



6 Z = arccos(e z ■ 5/|S|), 



(35) 



where e z is the nominal direction of the beam. 

Figure 7 plots the corresponding 5 Z versus A!T SC for each de- 
tector, and the numerical values are tabulated in Table 6. There 



LFI18M 5.94 7.19 

LFI18S 7.35 9.10 

LFI19M 6.83 8.49 

LFI19S 7.92 10.09 

LFI20M 9.43 11.56 

LFI20S 10.10 12.65 

LFI21M 9.46 11.59 

LFI21S 9.78 12.31 

LFI22M 6.96 8.62 

LFI22S 8.20 10.38 

LFI23M 5.68 6.90 

LFI23S 7.23 8.98 

70 GHz (mean) 7.91 9.82 

LFI24M 2.37 2.95 

LFI24S 1.43 1.98 

LFI25M 1.31 1.95 

LFI25S 0.88 1.49 

LFI26M 1.31 1.91 

LFI26S 0.90 1.44 

44 GHz (mean) 1.37 1.95 

LFI27M 8.54 11.11 

LFI27S 7.92 10.78 

LFI28M 8.66 11.25 

LFI28S 7.92 10.77 

30 GHz (mean) 8.26 10.98 

" Deflection in the polar axis of the dipole, in arcmin (Eq. 34). 
* Correction to the amplitude of the dipole, in ^/Kcmb (Eq. 35). 



is an evident linear correlation between the two quantities, with 
the best fit being 



1.2368 A7\ c + 0.2329. 



(36) 



The deflection angles are between 1.3' and 12.6', so that the 
largest correction is 13 fiK. It is interesting to note the slight 
level of asymmetry between the side and main arms of each 
feed-horn. This might lead to differences in corrections nearly 
as large as 1 fiK. An example is given by the two arms of LFI28. 
In contrast, the two arms of LFI21 show no asymmetries. 

As we already stated in Sect. 2.2, for this data release we 
have neglected the A!T SC (thus setting ||5|| = 1), but we have con- 
sidered the tilt effect given by 5 Z . The tilt was considered not 
only when using Eq. 25 (i.e., fitting timestreams with the ex- 
pected dipole signal), but also when we cleaned the calibrated 
data of the dipole itself (see point 7 in Sect. 3.1). 

4.3. Imperfect behaviour of the ADC 

As stated in Sect. 3, the calibration process also includes the 
conversion of the analogue detector voltages by the ADCs. This 
should be absorbed into the final radiometer calibration, but the 
ADCs can suffer a differential linearity problem, where the volt- 
age step between binary levels varies slightly. This can lead to 
a small but sharp deviation in the response curve as shown in 
Fig. 10 of Planck Collaboration III (2013). Since the slope of 
the response curve determines the effective calibration, a large 
gradient change can happen over a small range of voltage, po- 
tentially becoming the dominant source of gain uncertainty. 

Such an effect is seen in the LFI radiometers for some chan- 
nels, especially those with low back-end gain, where the quanti- 
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Fig. 7. Regression between dipole deflection (arcmin) and 
dipole amplitude correction ( /iK) for each radiometer. Numbers 
are taken from Table 6. Side arms are in black, main arms are in 
grey and have been shifted by 5' to improve readability. 



zation levels will be more evident. For this reason the radiome- 
ters of 44 GHz are the most affected, as shown in the upper panel 
of Fig. 8 for the worst affected LFI25S-1 1 channel (for more in- 
formation, refer to Planck Collaboration III 2013). 

The correction for this effect requires the reconstruction of 
the response curves, which can be inferred from the related dif- 
ference in behaviour between white noise and detector voltages. 
The reason is that the white noise, like the dipole and sky tem- 
perature fluctuations, is sensitive to the gradient of the response 
curve, whereas the detector voltages are just determined by the 
point on the curve. This is explained in Appendix A of Planck 
Collaboration III (2013) and an example for LFI25S-1 1 is shown 
in that paper together with estimates of ADC residual errors in 
terms of systematics at the map level. The implementation in the 
pipeline uses splines 11 , fitted to pre-calculated template curves 
of the non-linearities. The result of correcting the TOI data is 
shown in the lower panel of Fig. 8, where the dipole gains and 
inverse sky voltage and white noise are much more consistent. 

In the case where calibration used the level of the 4 K load, 
the uncorrected data would produce uncertainties directly in 
the calibration of the order of 5 % leading to stripes of about 
150 yuK, due to the solar dipole being incorrectly removed. 
Using a dipole based calibration would correct the ADC non- 
linearity, but the need to low-pass filter it for an acceptable 
signal-to-noise ratio would mean that gain variations shorter 
than the filter timescale were not corrected. This in turn would 
result in fine striping at the level of 20-30 fiK. Further, since 
the ADC non-linearity can differ between the two arms of a ra- 
diometer, there is also a direct impact on polarization through 
differential calibration errors. This is illustrated in Fig. 9, where 
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1 1 We use the spline functions provided by the GNU Scientific Library 
(Galassi et al. 2003). 



Fig. 8. The top panel shows the comparison for ADC uncor- 
rected LFI25S-1 1 diode data between the raw dipole gains (grey 
dots) and the reciprocal of total-power sky detector voltage 
(black line), the latter scaled to have the same mean level as the 
gains. Two highly significant departures can be seen around day 
160 and 560 together with lower level variations over most of the 
time range. The thick dark grey line shows inverse white noise 
estimates, again scaled to the same gain level, but filtered with a 
100 ring moving median filter to reduce the scatter. As expected, 
this follows the dipole gains. The lower panel shows the same 
thing, but after the ADC correction has been applied to the TOI 
data. 



we show the difference between maps made from the main and 
side arms of radiometer LFI25. 

It should be pointed out that the LFI and HFI ADC non- 
linearity problems are completely different, due to the differ- 
ent implementations of fast switching schemes. LFI use phase 
switching to swap sky and reference load signals, whereas HFI 
use a square wave modulation and AC coupling to avoid l/f 
contributions from the back-end electronics. More details about 
the ADCs used by HFI are found in Planck Collaboration VIII 
(2013). 



5. Accuracy of the calibration 

In this section we present results of consistency checks sensitive 
to calibration inaccuracies, as well as simulations which quan- 
tify the estimated level of calibration systematics in the data. 
Table 7 shows an overview of the elements that have been taken 
into account in estimating the calibration accuracy for LFI maps. 
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Table 7. Elements considered in deriving the accuracy of the calibration. 



Type of uncertainty 


Applies to 


Method used to assess the accuracy 


Reference 


Absolute 








Standard 


All sky 


Propagation of WMAP errors 


Sect. 5.1.1 


Zero level 


All sky 


Comparison with WMAP values 


Sect. 3.4 


Beam uncertainty 


All sky 


GRASP model of the beams 


Sect. 5.2.2 


Sidelobe convolution effect 


All sky 


Simulations 


Sect. 2.2 


Colour corrections 


Galactic areas 


Comparison of ground/flight bandpass leakages 


Sect. 2.3 


Relative 








Statistical/algorithmical errors . . . 


All sky 


Simulations 


Sect. 5.1 


Known systematics 


All sky 


Simulations 


Sect. 5.2 


Unknown systematics 


All sky 


Null tests 


Sect. 5.3 



Table 8. Accuracy in the calibration of LFI data. 



Type of uncertainty 
Absolute 

Standard" 

Zero level [^Kcmb] 

Beam uncertainty 

Sidelobe convolution effect 

Colour corrections 

Relative 

Statistical/algorithmical errors 6 [yuK CMB pixel -1 ] 

Known systematics'' 

Unknown systematics^ [/jKcmb pixel -1 ] 

Unknown systematics c (50 < I < 250) 

Total 

CMB areas^ [ //K C mb pixel - 1 ] 

Galactic region^ [/jK C mb pixel -1 ] 

Sum of absolute and relative errors^ 



Applies to 



30 GHz 



44 GHz 



70 GHz 



All sky 
All sky 
All sky 
All sky 
Galactic areas 



All sky 
All sky 
CMB areas 
Galactic region 
All sky 



All sky 



2.23 



0.25 
-300.84 

0.5% 

0.2% 
\a-2\ 0.1% 



4.3 
0.1% 

< 8.8 

< 17.0 
0.2% 



< 8.5 

< 38.5 
0.82% 



0.25 % 
-22.83 ± 0.78 

0.1% 

0.2% 
\a - 2| 0.3 % 



4.7 
0.1 % 

< 5.2 

< 9.8 
0.2% 



< 7.1 

< 13.7 
0.55 % 



0.25 % 
-28.09 ± 0.64 

0.3% 

0.2% 
\a - 2| 0.2 % 



6.5 
0.1% 

< 9.5 

< 13.1 
0.1% 



< 8.2 

< 16.8 
0.62 % 



" Error on the estimation of the calibration constant K iip (Eq. 27). 

* Peak-to-peak differences, as reported in Planck Collaboration III (2013). 

c Error in the estimation of the calibration constant K. 

d Scaled r.m.s. of the value of pixels in odd-odd survey half-difference maps minus the r.m.s. due to statistical noise. We used here JV side = 128 for 

30 and 44 GHz, and N Mc = 256 for 70 GHz. 
c Estimated from inter-channel comparisons on the cross-spectra. 

f Scaled r.m.s. of the value of pixels in odd-even survey half-difference maps minus the r.m.s. due to statistical noise. As above, we used here 

JV side = 1 28 for 30 and 44 GHz, and N siic = 256 for 70 GHz. 
s Sum of the error from the standard uncertainty and the square root of the squared sum of the following errors: (1) beam uncertainty; (2) sidelobe 

convolution effect; and (3) unknown systematics (as measured from the power spectrum at 50 < t < 250). 



Table 8 quantifies the impact of each element in the overall ac- 
curacy budget for the calibration. 



5.1. Statistical uncertainty 

5.1.1. Propagation of statistical errors from WMAP's dipole 
estimate 

As stated in Sect. 2.2, we have calibrated LFI data using the 
CMB dipole as measured by WMAP (Jarosik et al. 201 1). In our 
calibration we used both the speed of the Solar System in the 
CMB rest frame (vs un in Table 3) as well as the direction of mo- 
tion (/, b). In this section we estimate how the quoted errors on 
these parameters propagate through the calculation of the cali- 
bration constant K& v in Eq. 27. 



We ran the LFI calibration pipeline a number of times for a 
subset 12 of radiometers; for each run we changed either vs un or 
/ from the nominal value listed in Table 3, and we compared the 
new average level (^ ip ) of the calibration constant with respect 

to the level (^dip) of the nominal calibration. We also ran these 
tests without applying the iterative calibration code that cures 
the bias induced by the CMB and Galactic signal in dipole fitting 
(Mademoiselle; see point 5 in Sect. 3.1), in order to quantify 
the improvement in the stability of the calibration achieved by 
the iterative code. 

Figure 10 shows the results of our simulations in the case of 
LFI27M (results for the other radiometers are in nearly perfect 



12 Specifically, we considered LFI18M, LFI18S (70GHz), LFI24M, 
LFI24S (44 GHz), LFI27M, and LFI28M (30 GHz). 
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Fig. 9. Three degree smoothed map of the difference between the 
main and side arm of LFI25 before (top) and after (bottom) ADC 
non-linearity correction. Stripes can seen which are stronger in 
the region of the solar dipole due to residuals from calibration 
uncertainties. This relative uncertainty between radiometer arms 
directly affects polarization as shown here, since the main and 
side arm difference represents a pseudo-polarization. 



agreement with this). Mademoiselle's ability to reduce the dis- 
crepancy (^ip) _ (^ di p) ' s cleax if one compares the error on K^ v 
with and without using Mademoiselle. Considering WMAP's 
uncertainties on vs U n and I and the results of our simulations, we 
can estimate that the relative error cr^/K^ due to the uncer- 
tainties in vsun, I, and b is 

— ^ = 0.19 % (vsun) + 0.05 % (/) + 0.01 % (b) 

K,:v, (37) 

= 0.25 %, 

where the error for b has been estimated by simple proportional- 
ity using the error on I. (These numbers are the same for all the 
radiometers we analysed.) 

5.1 .2. Simulation of noise in dipole fitting 

We have run a set of simulations to assess how much the statis- 
tical noise in the radiometric signal affects gain reconstruction 
using the CMB dipole signal. We first decalibrate the input map 
or dipole with a fiducial gain, the 4 K calibration, and then we 
simulate how the dipole calibration pipeline reconstructs those 
gains, i.e., we add the statistical noise, cut the low-dipole regions 
and then apply a smoothing filter. 

We simulated separately the impact of the statistical noise 
acting on the modelled dipole and the Galaxy/CMB contribu- 
tion; we model the input dipole using the WMAP direction and 
amplitude for the cosmological dipole (see Table 3) and NASA 
JPL Horizons Planck velocity for the orbital dipole. Moreover, 




Error on the Solar System speed [km/s] 




n 1 1 — — i 1 1 r 

-1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 



Error on the dipole axis angle [°] 

Fig. 10. Variation in the average level of the calibration con- 
stants of LFI27M, with (black line) or without (grey line) ap- 
plying the Mademoiselle iterative code, as a function of the 
error in the Solar System speed (top) and the displacement of 
the dipole axis (bottom). The two plots are normalized so that 
all the lines go through the point (0,0). The width of the verti- 
cal grey band corresponds to the error quoted by Jarosik et al. 
(201 1) for vs un (top) and / (bottom). 



we used Planck full frequency maps as the input Galactic and 
dipole-removed CMB signal. 

We use a fiducial gain in order to make our simulation 
as realistic as possible, including the impact of the calibration 
pipeline on the underlying "true" gain. The 4K calibration (see 
Sect. 3.2) offers a good choice, because it also includes realis- 
tic fluctuations on shorter timescales, which are expected to be 
smoothed out by the dipole calibration algorithm. 
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Fig. 11. Power spectrum of the simulated artefacts caused by 
errors in gain reconstruction (see Sect. 5.1), compared with 
Planck's, best-fit TT power spectrum model. 



we produced time-ordered data by summing Galactic, dipole, 
and CMB signals, as well as 1// noise. To do this we considered 
Planck's scanning strategy and the measured characteristics of 
each LFI's detector. We de-calibrated the TODs using the same 
gains used to calibrate flight data. We then produced two sets of 
maps: 

1. the first set uses these TODs and produces maps, which are 
immune from the ADC effects; 

2. the second set filters the TODs using the inverse of the 
response function R in order to simulate the ADC non- 
linearity, and then uses the same procedure used for flight 
data to estimate a new response function, R' (obviously, 
R * R'); finally the ADC correction code is applied to the 
data, and new maps are produced. 

We considered the difference between the two sets of maps to be 
a reasonable estimate of the level of errors due to the uncertain- 
ties in the procedure for estimating ADC non-linearities. These 
differences reveal that such uncertainties introduce an error in 
the estimation of the calibration constant K (Eq. 1) which is of 
the order of 0.1%. 



After decalibrating the input signal with the fiducial gain, 
we add a realization of the expected statistical gain error related 
to the dipole fit due to white noise and to Galaxy masking. We 
model the statistical gain error according to the following equa- 
tion: 

o~ c ,i _ V2o-q 



1 + 



A. 

,/kne 



(38) 



where i is an index identifying a pointing period, cr c ,/G, is the 
relative error on the gain G-„ cr is the noise-equivalent temper- 
ature of the radiometer, D, is the amplitude of the dipole signal 
during pointing period i (in Rayleigh-Jeans units), r, is the in- 
tegration time of pointing period i, < /7mask,t < 1 is a dimen- 
sionless coefficient that takes into account the loss of integra- 
tion time due to Galaxy masking within pointing period i, /o is 
the frequency of the dipole sinusoid seen by Planck's scanning 
strategy (/o = 16.7 mHz), and f^ nee and a are the constants that 
characterize the 1 // noise for a given radiometer. The output of 
this step is equivalent to the "raw" dipole calibration as com- 
puted by the Mademoiselle iterative calibration software (see 
Sect. 3.1, point 5, and the Planck Collaboration II 2013). 

The last step is to cut the low-dipole regions, replace them 
with a straight line, and then heavily smooth the data with a 
moving average filter of about 200 pointing periods. This is a 
simplified version of the wavelet-based filtering in the pipeline 
(Sect. 3.1, point 6). We run destriping on the output timelines 
and produce surveys and nominal mission maps of the noise for 
each frequency. Figure 1 1 shows the power spectrum of the arte- 
facts produced by this simulation for the same timespan covered 
by this Planck data release. 

5.2. Impact of known systematics 
5.2.1. ADC non-linearities 

We have run a set of simulations on the non-linearity of the 
ADCs (see Sect. 4.3) in order to quantify its impact on the 
calibration of LFI. The simulations are described by Planck 
Collaboration III (2013), but we briefly recall the procedure here. 

After we characterized the non-linearity of the ADCs to esti- 
mate the response curve /?( V/) of the ADC for the 44 converters, 



5.2.2. Beam uncertainties 

In Table 8 we quantify the error due to the fact that we con- 
sidered monochromatic beams when estimating the vector S in 
Eq. A. 5. We have estimated the beam response of LFI's radiome- 
ters (see Sect. 4.2) over their bandwidth using the GRASP soft- 
ware (see Planck Collaboration III 2013, for more information 
about this) and we have computed a new set of vectors S'(v), 
which of course depend on the frequency v. The errors we report 
are equal to the average of the value ||S|| - [|<S'0)>|| for the 22 
radiometers, grouped according to their centre frequency. 

5.3. Impact of unknown systematics 

To estimate the impact of all those systematics which affect the 
calibration, but for which we do not have simulations or direct 
measurements at hand, we relied on internal consistency checks, 
namely inter-channel comparisons and null tests. We explain the 
logic behind each analysis as well as the results in the next para- 
graphs. 

5.3.1. Inter-channel calibration consistency 

A sensitive way to test the relative calibration among the differ- 
ent LFI radiometers is to make independent maps and compare 
their angular power spectra in the vicinity of the CMB acous- 
tic peak. This has the advantage of testing the calibration on a 
source with the relevant amplitude and spectrum for our cosmol- 
ogy goals. To reduce the effects of noise bias on the power spec- 
trum estimation, we intercompare the radiometers using cross- 
spectra; since the noise is quite uncorrected among the radiome- 
ters, we obtain clean spectra this way. In order to compare any 
pair of radiometers, we first mask the Galaxy, point sources, and 
all unobserved pixels. We then calculate the cross-spectrum of 
maps made by each radiometer with maps of a third radiometer, 
and take the ratio. By comparing all possible combinations we 
average over any small impact from the choice of the third map. 
For example: 



C 



LFIlSSxother 
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^LFI19Mxother ' 



(39) 
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Fig. 12. Intercomparison of the relative calibration of the LFI ra- 
diometers, as measured on the CMB acoustic peak by averaging 
the power spectrum from 100 < I < 250. Note that we display 
both main arm (slightly to the left of each dotted, vertical line) 
and side arm (to the right) radiometers for each horn, as they are 
all calibrated independently. We also include the full frequency 
maps, labelled as /30, /44, and /70. LFI 70 GHz full frequency 
at 70 % sky coverage is taken as reference. 



In this case, "other" runs over LFI18M, LFI18S, and both ra- 
diometers of horns 20, 21, 22, and 23. 

We show the results of this in Fig. 12, which captures the rel- 
ative gain deviation among all the radiometers, based on the first 
acoustic peak. The horizontal axis is the horn number, with main 
and side arm radiometers spaced before and after the horn num- 
ber. Figure 12 also includes the gain deviations averaged over the 
full frequency, labelled as /30, /44 and /70. Interchannel con- 
sistency is everywhere within 1 % of the mean at 44 and 70 GHz, 
while 30 GHz is severely affected by residual Galactic signal, as 
highlighted by the strong reduction between the 60 % and 70 % 
sky coverage masks. This analysis includes corrections for the 
beams, the scanning strategy and mask deconvolution. 

Figure 13 shows the results of a similar comparison between 
surveys, in this case comparing the full frequency estimates of 
the acoustic peak for surveys 1, 2, and 3. The full frequency 
comparison is even more striking, with each band remaining 
within ±0.15 % over the three independent surveys. Of course, 
many systematics which plague calibration may be common to 
all channels, particularly within a frequency cohort; however, 
1 // noise, receiver degradation, and certain other sources of sys- 
tematic uncertainties in the temperature calibration seem to be 
controlled to below 0.15 % . 

5.3.2. Null tests 

We have developed a number of tests to highlight possible issues 
related to instrumental systematic effects which are not properly 
corrected within the pipeline, or are related to known changes 
in the operating conditions of the instrument or to intrinsic in- 
strument properties coupled with the sky, like stray-light from 
sidelobes. Such tests are also discussed by Planck Collaboration 
III (2013); here we concentrate on how those tests can be used 
to estimate the impact of calibration errors in LFI maps. 




Survey number 

Fig. 13. LFI cross-spectra between surveys. As in Fig. 12, ratios 
of the power spectra are computed in the range 50 < t < 250. 



The idea behind the many tests we ran for this Planck data 
release is to compare pairs of distinct estimates of the sky signal, 
that should ideally be the same, such as: 

- comparing the sky signal measured during two surveys (not 
necessarily consecutive), we call this kind of test survey- 
difference null tests; 

- comparing the sky signal measured by two radiometers, we 
call this a radiometer-difference null test; 

- comparing the map produced using only data taken in the 
first half of each pointing period with data taken in the other 
half. This is a half-ring null test. 

Clearly, it is not trivial to determine the cause of each discrep- 
ancy in each comparison. The following criteria are useful in 
interpreting the results of these null tests: 

1. The results of radiometer-difference null tests can reveal 
artefacts due to imbalances in the frequency response, as 
well as to their different optical systematics. Such optical ef- 
fects can be caused by an error in pointing reconstruction 
(which is therefore unrelated to calibration), but also by two 
different ways of picking up the astrophysical signal (includ- 
ing the dipole) through the sidelobes. The latter effect can 
be critical for calibration, as the subtraction of an improper 
dipole signal leads to unwanted residuals in the map. 

2. Survey-difference null tests are not affected by frequency re- 
sponse mismatches. However, they can be affected by the op- 
tical effects described above. Because of the scanning strat- 
egy employed by the Planck spacecraft, such systematics 
will be present only if an odd number survey is compared 
with an even survey, or vice versa, since odd surveys share 
the same scanning strategy, and the same applies for even 
surveys. (This is true for the first four surveys, but not for 
subsequent surveys, where the precession phase angle was 
changed.) 

3. Half-ring null tests are able to spot systematics at high fre- 
quency (more than 1 /60 Hz, as this is the spinning frequency 
of the Planck spacecraft). Thus, calibration errors are un- 
likely to affect this kind of test, since we expect the gain of 
the radiometers to change on much longer time scales (days 
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Fig. 14. Survey difference maps at 30 GHz. Top: difference be- 
tween survey 1 and 2. Middle: difference between survey 1 and 
3. Note that the 1-2 difference map has sharper features, due 
to the different scanning strategy between odd and even sur- 
veys that prevents optical systematics from being cancelled in 
the subtraction. Bottom: the mask used to distinguish between 
CMB-dominated regions and Galactic regions when computing 
the numbers reported in Table 8 (30 GHz case). 

or weeks). In fact, we expect half -ring null test to be sensitive 
only to the level of unavoidable white noise in the data. 

For all the reasons stated above, we believe that a fair up- 
per limit of the calibration errors that are unrelated to optical ef- 
fects can be deduced by studying odd-odd and even-even survey- 
difference null tests, as they are not affected by pointing errors, 
imbalance in frequency response, or "fast" systematics. Even- 
odd and odd-even differences are affected by systematics caused 
by optics as well. Therefore, these should be larger than the for- 
mer and can also be used to produce upper limits for the calibra- 
tion accuracy. Figure 14 shows two examples. 

We have used odd-odd survey difference maps to produce 
the numbers reported in the row "Unknown systematics" in 
Table 8. These numbers are the r.m.s. value of the pixels, after 
the maps have been degraded to a resolution which roughly co- 
incides with the FWHM of the LFI beams. We subtracted from 
these values the part due to statistical noise, estimated from half- 



ring difference maps using the procedure described by Planck 
Collaboration III (2013). As said above, these r.m.s. values rep- 
resent upper limits for calibration inaccuracy, because some of 
the systematics affecting survey difference maps are unrelated to 
calibration. Values under the "Total" group have been calculated 
from odd-even survey difference maps; as noted above, they in- 
clude optical systematics as well. 

We show the power spectra of the even-odd and odd-odd null 
tests in Fig. 15, where we compare them with the estimated sta- 
tistical errors of the gains (Sect. 5.1) and the spectra of half-ring 
maps. The latter provide an estimate of the level of statistical 
noise in the maps and generally agree with survey difference 
spectra for t > 20. This means that for a wide range of multi- 
poles we do not expect any significant contribution of calibration 
systematics on the maps. For i < 20 there is some disagreement 
for the 30 GHz maps. This might indicate the presence of resid- 
ual optical systematics in the data. In fact, the chief reason that 
induced us to use a different calibration for 30 GHz channels (us- 
ing the 4 K total -power output, see Sect. 3.2) was to minimize the 
impact of optical systematics on the calibration of these channels 
(which suffer from a relatively high level of sidelobes). 

5.4. Checks with external calibrations 

5.4.1 . Measure of point source flux densities 

In order to check the accuracy of the calibration between the 
different RCAs, and to test the colour corrections and central 
frequencies of the bands, we look at three bright radio sources 
in the LFI maps: Orion A (M42); Taurus A; and Cassiopeia A 
(none of which are expected to vary). The emission in these three 
sources is due to synchrotron or free-free emission, such that 
they have very different spectra from the CMB. This means that 
colour corrections are important at the level of about 3-10% 
depending on the RCA, which is in excess of the uncertainties 
on the flux densities of these sources. 

Comparisons between ERCSC flux densities (Planck 
Collaboration VII 2011) at 30, 44, and 70 GHz and inter- 
polated or extrapolated flux densities from Karl G. Jansky 
Very Large Array (VLA) follow-up observations have pre- 
viously shown good statistical consistency at the 3 % level 
(Planck Collaboration XIV 2011; Kurinsky et al. 2012; Planck 
Collaboration XXVIII 2013). A comparison has also been made 
between the Plcmck-AICA Coeval Observations (PACO) and 
PCCS catalogues in Planck Collaboration XXVIH (2013), which 
finds consistency at the level of around 4 %. However, we ide- 
ally want to test the calibration much more accurately than this, 
which can be done internally and with comparison to the WMAP 
results. Additionally, the sources that we examine here are sig- 
nificantly less variable than the typical source in the VLA and 
PACO tests. 

We use the Planck single detector maps for each survey. We 
first average the main and side maps to create individual RCA 
maps per survey, which has the effect of cancelling out the po- 
larization such that the maps are of total intensity only. We look 
at the individual survey maps, and the average of surveys 1-3. 
We also make use of the WMAP-9 deconvolved Stokes / sky 
maps (Bennett et al. 2012). 

We measure source emission from the maps using aperture 
photometry, as described in Planck Collaboration XX (2011), 
with an aperture of 70' radius and a background annulus of 80'- 
100'. By keeping the aperture and annuli sizes constant rather 
than changing them to reflect the resolution of the instrument, 
we ensure that we are measuring the same source extent and 
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Fig. 15. Spectrum of the half-difference between Sky Survey (SS) 1 and 2 ((SSI - SS2)/2, odd-even, frame A) and (SSI - SS3)/2 
(odd-odd, frame B), compared with the half ring spectra for the same time period at the three centre frequencies of LFI and the 
simulated errors due to gain reconstruction (Sect. 5.1). Spectra have been binned into five samples each, in order to remove high- 
frequency noise and highlight the general trend. At high multipoles (I > 20) there is a good match between survey differences 
and half-rings, which means that discrepancies between survey differences are mainly due to statistical noise. Note that the biggest 
discrepancy between survey differences and halfrings happens for 30 GHz channels when comparing SSI with SS2 (odd-even); 
this suggests that optical effects (sidelobes) are causing the discrepancy. Gain errors are calculated over the whole mission and are 
negligible except for the lowest multipoles of 44 GHz and 70 GHz spectra, where they probably account for the large-scale residuals 
found in survey difference maps. 



subtracting the same background emission at all frequencies, 
which may vary with distance due to the structure of the ex- 
tended source, or background emission surrounding it. A cor- 
rection factor for the source flux that falls outside of the aperture 
is applied as per Planck Collaboration XXVIII (2013); this fac- 
tor is non-negligible for the WMAP data points. The measured 
emission, in 7cmb, is converted to Trj, and then into Jansky at 
the nominal frequency (28.4, 44.1 or 70.4 GHz as appropriate). 
The uncertainty in the flux densities is calculated using the r.m.s. 
of the pixels in the background annulus. 

The properties of the sources are summarised in Table 9. The 
flux densities for the sources are the aperture photometry flux 
densities from PCCS; however as the PCCS only contains the 
flux densities for the RCA-averaged maps these are only used 
here for consistency checks. We note that there is a difference in 
reported flux densities from WMAP between Page et al. (2007) 
and Weiland et al. (201 1) of up to around 8 %, due to the method 
used: Page et al. (2007) used aperture photometry, and Weiland 
et al. (201 1) fitted the beams to the sources. As such, care needs 
to be taken when comparing the numbers from the different anal- 
yses. 

In order to calculate the necessary colour corrections, we it- 
eratively calculate the amplitude (at 33 GHz) and spectral index 
of the source based on the Planck and WMAP aperture photom- 



etry flux densities, the values for which are quoted in Table 9. 
We then use this spectral index to colour-correct the flux densi- 
ties per RCA using quadratic fits to the colour correction values 
given in Sect. 2.3. The average flux densities, and the standard 
deviation calculated between the RCAs, both before and after 
colour correction, are given in Table 9. The percentage of the 
standard deviation to the mean is also given. When the colour 
corrections are applied, the standard deviations significantly de- 
crease, giving confidence that the colour corrections are improv- 
ing the consistency between the RCAs. 

Figure 16 shows the spectra of the three sources (left-hand 
column) and the flux densities from individual survey maps for 
the different RCAs (right-hand column). For Orion A, two sets of 
points differ significantly from the spectrum from aperture pho- 
tometry. At WMAP 93.5 GHz, the measured flux density is much 
higher as it includes dust emission, which is not present at the 
LFI frequencies and as such is not considered here. Additionally, 
the PCCS aperflux is much steeper (a = -0.21) than the best 
fit to the aperture photometry here; this is most likely due to 
the different sized apertures used, which depend on frequency 
for PCCS, meaning that it will systematically exclude the more 
diffuse emission from the source that is included in this anal- 
ysis. This is particularly relevant for Orion A due to its more 
complex morphology. For Cas A, we find very good agreement 
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Fig. 16. Spectra (left) and per-survey and per-RCA measurements (right) for Orion A (top), Cas A (middle) and Tau A (bottom). 
For the spectra, the flux densities measured using aperture photometry are shown for LFI and WMAP, as well as the best-fitting 
power-law to the measurements (solid line) and a ±1 % range (dotted line). The fit and data points from Weiland et al. (2011) are 
shown for Cas A and Orion A, and the Planck Catalogue of Compact Sources (PCCS) flux densities for all three sources are shown. 
For the per-RCA measurements, RCAs 18-23 are at 70 GHz, 24-26 at 44 GHz and 27-28 at 30 GHz. For each RCA, the average 
from the three surveys is shown, followed by the measurements from each survey. 



between the results. For Tau A, we find a slightly flatter spec- 
trum (-0.287 + 0.013) than the PCCS and Weiland et al. (201 1) 
(-0.302 ± 0.005), and other analyses also report steeper spectra, 
e.g. Hafez et al. (2008) find -0.32 + 0.01 at 30 GHz, and Maci'as- 
Perez et al. (2010) find -0.296 + 0.006; this is again likely due 
to the different aperture sizes considered here. 

Although aperture size affects the comparison of the absolute 
flux densities, it does not affect the results between the different 
Planck RCAs at the same frequencies, or between the different 
survey maps for the same RCAs. There is general consistency 
for all RCAs between the surveys, within the uncertainties. Note 



that some of the sources were not observed in every survey due 
to the survey definitions, for example Cas A was not observed by 
RCAs 25 and 26 in survey 2. However, we do find some system- 
atic differences in the results from different RCAs. At 28.4 GHz, 
RCA 27 is consistently lower than RCA 28 at the level of 2 %; 
this difference does not change when colour corrections are ap- 
plied as the corrections for those two RCAs are so similar. The 
differences at 44 GHz are 1 .5 %, with RCA 24 systematically 
slightly higher than the 25 and 26. This might be due to the 
beams, as RCA 24 has a smaller beam size than RCAs 25 and 
26, however this difference is expected to be smaller than the 
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Table 9. Properties of three compact calibrators used to check 
the calibration: Orion A (M42); Cas A; and Tau A (Ml) 



Property 


Orion A 


Cas A 


Tau A 


Position (/, b) 


(209?01, -19T38) (111?74, -02! 14) (184?56, -5?78) 


c 28.4 GHz rT,,na 

j pccs LJyi 

r 44.1 GHz r, , 

J pccs L-<yj 

s™r [Jy] 


365 + 2 
341 ±3 
301 ±3 


192.0 ± 1.4 
137.0 ± 1.2 
102.3 ± 0.5 


352.4 ± 0.9 
301.8+1.3 
260.0 ± 0.9 


A" 

a" 


370.5 ± 1.7 
-0.026 ±0.014 


179.9 + 0.8 
-0.668 ± 0.014 


319.1 ± 1.0 
-0.287 ± 0.013 


sl™ GHz UyY 
s^ GHz [Jy] 
s™ GHz [Jy] 


384 ± 4(1.0%) 
373 ± 2(0.5%) 
372 ± 10 (2.7 %) 


204 ±3 (1.5%) 
151 + 2(1.2%) 
110 + 4(3.9%) 


340+ 4(1.3%) 
299+ 2(0.5%) 
255 + 11 (4.2%) 


C 28.4 GHz rT„l<' 
°aper,cc l J }i 
r44.1GHz ry , 
u aper,cc V J J i 
r 70.4 GHz r| , 
a aper.cc L-VJ 


375 ± 4(1.0%) 
370 ± 3(0.7%) 
364 ± 7(2.0%) 


198 + 3 (1.5%) 
149 + 2(1.5%) 
106 ± 3 (2.8 %) 


334 + 4 (1.2%) 
296+ 1.2 (0.4%) 
249 + 4 (1.8%) 



" Flux densities from PCCS, without colour correction. 

* Coefficients for the fitted power law, 5 = A(y/33 GHz) ff . 

c Average and standard deviation of the flux densities measured by 
aperture photometry. 

d Average and standard deviation of the colour-corrected aperture pho- 
tometry flux densities. 



Table 10. Average ratio of flux densities for RCA pairs. 



Ratio 


6C a= -i 


8C„ = o 


Corrected 
Discrepancy discrepancy" 


70 GHz 










S RCA19 1 'S RCA18 • • 


. 0.073 


0.049 


-1 % 


6% 


5rcA2o/5'rca18 • • 


. 0.047 


0.032 


0% 


4% 


S'rcA2i/S , RCA18 • ■ 


. 0.025 


0.017 


0% 


2% 


5rcA22/5rca18 • ■ 


. -0.055 


-0.035 


7% 


1% 


5rcA23/5'rca18 • • 


. -0.024 


-0.013 


2% 


0% 


44 GHz 










S , RCA24/S'rCA26 • • 


. -0.015 


-0.011 


0% 


1% 


5 , RCA25/5rca26 • • 


. -0.007 


-0.005 


1% 


1% 


30 GHz 










5 , RCA27/5rca28 • • 


. 0.001 


0.000 


2% 


2% 



Value of the column "Discrepancy" minus the difference of colour 
corrections for a = -1 (see Table 1). 



effect seen here. At 70 GHz, there are indications that RCA 18 
is consistently high and RCA 19 is consistently low, with differ- 
ences in the flux density between those two RCAs of up to 6 % 
(for Orion A). 

In order to better characterize these differences, we have 
studied how much the ratio Sa/Sb, the average ratio between 
compact source flux densities calculated using RCAs A and B, 
differs from unity. In addition to aperture photometry, we have 
estimated this ratio using two other methods: 

1 . We used the template fitting code from Davies et al. (2006), 
run in a simplified way (i.e., without making use of a covari- 
ance matrix for noise or the CMB, such that it essentially 
does a least-squares fit), to compare "template" maps of 



RCA 28 (30 GHz), RCA 26 (44 GHz) and RCA 18 (70 GHz) 
with the other maps at the same frequencies. The maps have 
been smoothed to 1°, and a constant offset is used as a second 
template to remove any differences in the offsets between the 
template and the map. This provides an estimate of the ratio 
Sa/Sb for the all-sky coefficients between RCAs. We have 
run this analysis on the full, unmasked sky and the Galactic 
plane (\b\ < 10°), which yields consistent results with those 
obtained through aperture photometry; in addition we have 
used the WMAP KQ75 mask on 3°-smoothed maps and find 
consistent results at 30 GHz, but better agreement at 44 and 
70 GHz, where the maps are more dominated by the CMB. 
2. We combined pairs of RCA maps according to the formula 
niA - yniB, with twa and thb the arrays of pixels for the two 
maps, and visually compared the three sources considered 
above. We found those values of y that either made asym- 
metries in the residual beam patterns around point sources 
disappear (30, 44 GHz) or the source difference is within 
the noise uncertainty (70 GHz). We then assumed that y « 
Sa/Sb- Although this method is rather qualitative, the val- 
ues we obtained differ from aperture photometry estimates 
by only about 1 %. 

Table 10 shows a comparison of the average ratio between 
source densities calculated using two RCAs, e.g., RCA 18 vs. 
RCA 19. The table reports an estimate of the percentage discrep- 
ancy between the channels (calculated simply by averaging the 
discrepancies found by aperture photometry, template fitting and 
visual inspection of sources), as well as the discrepancy minus 
the expected difference in colour corrections between the RCAs 
(assuming a = -1; the difference of colour corrections for C a= o 
are also shown in the table for comparison). The expected un- 
certainties on the colour corrections for an RCA for a = - 1 are 
0.42 %, 1 .38 % and 1 .53 % for 30, 44, and 70 GHz, respectively, 
such that the effects we are seeing here are 4cr at 30 and 70 GHz. 

As the consistency between the RCAs is seen to be signif- 
icantly better than that measured here when comparing CMB 
emission between the RCAs (e.g., Sect. 5.3.1), we conclude that 
the effects seen here are likely due to differences in bandpasses 
or the intermediate beams between the RCAs, although neither 
of these is expected to be sufficiently large to explain the differ- 
ences seen here. These effects will be investigated further prior 
to the next release of Planck data. 

5.4.2. Planet flux densities 

Planets can be used as a source for calibration cross-checks. In 
particular, Jupiter has several advantages: its signal is up to fifty 
times higher than the receiver noise per scan; it is visible about 
twice a year; it has a relatively well known spectrum; and at 
LFI's resolution it is point-like. Table 11 lists the dates when 
planets crossed LFI beams, and Fig. 17 gives a visual timeline 
of these events. 

The model for the time dependent signal produced by a 
planet in a timeline (Cremonese et al. 2002; Maris et al. 2007) 
is: 

/p^AT 1 f dfl[B(x b )/p fl (*b-U(» i )Xpft-)) + /bgfe)], (40) 

where / Pi o is the brightness distribution across the projected 
planet disk, x v {t\) represents the direction in which a planet is 
seen from LFI at the sampling time r,-, U(f,) is the matrix that 
converts pointing from the Ecliptic reference frame to the beam 
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Fig. 17. Timeline of planet crossings. Each line shows the times when the planet fell within the beam of at least one of the 
radiometers at 30, 44, or 70 GHz. Gray bands at the edges mark the start of the survey and the end of the period considered in the 
2013 data release. Vertical grey lines mark the separation between sky surveys. 



Table 11. Visibility times of planets". 



limb darkening. Therefore, Eq. (40) reduces to: 





30 GHz 


44 GHz 


70 GHz 


Mars 


09/10/18-21 
10/04/17-20 


09/10/18-21 
09/10/29-11/02 
10/04/09-12 
10/04/18-20 


09/10/20-24 
10/04/15-18 


Jupiter 


09/10/31-11/02 
10/06/30-07/03 


09/10/24-27 
09/10/31-11/02 
10/06/30-07/02 
10/07/08-12 


09/10/29-11/01 
10/07/01-05 


Saturn 


09/12/31-10/01/03 
10/06/20-22 


09/12/31-10/01/03 
10/01/08-11 
10/06/10-13 
10/06/20-22 


10/01/02-05 
10/06/18-21 


Uranus 


09/12/12-15 
10/06/27-30 


09/12/13-15 
10/06/27-29 
10/07/06-08 


09/12/10-13 
10/06/29-07/02 


Neptune 


09/11/07 

10/05/15-17 

10/11/09-11 


09/10/31-11/02 

09/11/07-09 

10/05/15 

10/05/23-25 

10/11/10 


09/11/05-07 
10/05/16-19 
10/11/08-10 



Dates are in the format YY/MM/DD. Only the times between 2009 
Aug 13 and 2010 Nov 28 are reported here. 



reference frame, and 4g(f;) is the observed background. We as- 
sume that the radiometers have an elliptical Gaussian beam pat- 
tern, and we neglect complications such as differential temper- 
ature distributions across the planetary disk, phase effects, and 



y P,0 



A b 



(1 -/cos6» obs )exp 



+ /bg, (41) 



where y p = \Jx p , b is the beam FWHM, / is the planetary disk 
flattening, t9 bs is the observer planeto-centric latitude, and A is 
the distance between the planet and the observer. The diagonal 
matrix C = £> 2 /81og2diag(e, 1/e) represents the beam shape, 
whose ellipticity is e. The time dependency of I ? , 7 Pi o, A, # D bs, 
U, and Xp are omitted. The equation assumes a beam-averaged 
planet brightness 7 p o calculated for a fiducial planet-observer 
distance Ao and beam FWHM bo, in order to account for beam- 
to-beam variability. 

The model is used to fit 7 p o onto the timelines of each 
feed-horn, taken separately, by using a least-squares procedure. 
Samples for each timeline are accepted as valid data of the fit 
if their instantaneous angular distance from the planet is smaller 
than 5° at the epoch of observation. The background is assumed 
to be constant over the acceptance disc and is left as a free pa- 
rameter. In order to properly destripe the data, we have used the 
baselines calculated by Madam (Keihanen et al. 2010), the stan- 
dard map-making tool used to produce the LFI maps. 

The histogram of the data in the timelines can be approxi- 
mated as the sum of two Gaussian distributions, one due to the 
signal and one due to the background. The histogram of differ- 
ences of timelines and the fitted model overlap very well, one 
of the peaks assessing proper separation of the planetary sig- 
nal. The width of the residual peak is an estimate of the con- 
fusion noise, dominated by white noise and CMB fluctuations. 
Table 12 reports the measurements of the brightness tempera- 
ture for a number of planets. The two quoted errors are the max- 
imum intra-survey difference with respect to the median of the 
measures and the confusion noise measured on the fitted data. 

Jupiter is the only planet bright enough to be detected in the 
timeline without averaging. This is the reason why the random 
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Fig. 18. The Jupiter spectrum for LFI compared to the WMAP 
spectrum (Weiland et al. 2011). In both cases, measurement er- 
rors are smaller than the size of the symbols. 



Table 12. Averaged brightness temperatures for planets" . 







30 GHz 






44 GHz 






70 GHz 




Mars 


183 


± 1 ± 


4.1 


187 


± 10 ± 


4.1 


183 


± 3.5 ± 


2.3 


Jupiter 


134 


± 0.5 ± 


0.2 


148 


± 0.5 ± 


0.2 


157 


± 0.5 ± 


0.2 


Saturn 


121 


± 3.5 ± 


0.9 


128 


± 3.0 ± 


1.0 


131 


± 2.5 ± 


0.5 


Uranus 


190 


± 133 ± 


23 


230 


± 10.5 ± 


29 


138 


± 7.5 ± 


13 


Neptune 


79 


± 5 ± 


74 


74 


±22 ± 


76 


101 


±11 ± 


32 



Values are in Kelvin. The first error represents the intra-survey differ- 
ence, the second one the confusion noise. 



error is much smaller than the maximal survey-by-survey vari- 
ation. Such variation can be ascribed to a number of small sys- 
tematics not completely understood yet. One of them is the fact 
that planets have been extracted as soon as they were observed, 
so that different observations have had slightly different calibra- 
tion procedures. Despite having a slightly worse signal-to-noise 
ratio (S/N), the same is true for Saturn. However, we did not 
introduce any correction for the effect of Saturn's rings yet. In 
the same manner the observations of Mars are not corrected for 
diurnal variability due to the rotation of the planet. Uranus and 
Neptune have quite low S/N, in particular at 30 GHz, and there- 
fore they are very sensitive to the confusion noise, background, 
uncertainties in beam orientations, and aperture corrections. 

Figure 18 compares the Jupiter spectrum from LFI averaged 
through the surveys (diamonds) with the WMAP measurements 
(Weiland et al. 201 1). It is possible to se the excellent, albeit not 
perfect, agreement between the two spectra. However, we do not 
use such data to attempt a measurement of the quality of the LFI 
calibration: we feel that at this stage the uncertainties in the pub- 
lished data (seasonal variability of the source, inconsistencies in 
literature among measurements of the same data) and uncertain- 
ties in the analysis procedure would make such assessment too 
rough to be of any real use. 



6. Conclusions 

In this paper we have described the procedures used to calibrate 
Planck/LFFs data and the ways we used to assess the accuracy 
of our calibration. We have shown that the level of consistency 
in the maps is quite good, since discrepancies are roughly 0.6 % 
at 44 and 70 GHz, and 0.8 % at 30 GHz. 

However, we believe there are areas where we can signifi- 
cantly improve our understanding of the instrument and there- 
fore produce a better calibration, e.g., modelling of far sidelobes 
(see sects. 4.1 and 4.2) and proper modeling of variations in the 
noise temperature of the radiometers (Sect. 3.2). We therefore 
need to concentrate our efforts in improving our algorithms and 
analysis methods for the next Planck data delivery, in order also 
to provide a sound polarization calibration. 

An important point which is still missing is the production 
of a full characterization of the dipole that is independent of 
WMAP. This would in turn allow us to calibrate the Planck/LFl 
radiometers absolutely. As explained in Sect. 3.3, we were able 
to obtain some encouraging results, but not good enough to avoid 
using WMAP's characterization of the dipole in our calibration 
pipeline. We aim to do so for the next release. 

Another important point is to fully characterize the impact of 
optical systematics on the calibration. Our efficient dipole con- 
volver (Sect. 4.2) is a testimony to the work we have put in this 
task, but much more work is needed to fully understand the im- 
pact of far sidelobes on the calibration. We plan to run simula- 
tions to assess how much our uncertainty in the knowledge of 
the optical parameters of the beams impacts the estimation of 
gains and the production of calibrated timelines, and to further 
improve our model of the beams themselves. 

Finally, we are determined to understand the origin of all 
the "unknown systematics" listed in the corresponding row in 
Table 8, by means of simulations and more refined models of 
the radiometer. 
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Appendix A: Fast convolution of the beam with the 
dipole 

In this section we provide the details of the model used to con- 
volve the dipole signal used in the calibration with the beam re- 
sponse over the full sky (Sect. 4.2). 



A. 1 . Analytical model 

Our starting point is Eq. 25, which relates the dipole signal (as 
measured by a pencil beam) with the combined velocity of the 
spacecraft and of the Solar System with respect to the CMB rest 
frame. Following Peebles & Wilkinson (1968), we rewrite Eq. 25 
as 

= x(fd -De+Ux- D e ) 2 - l-\D E \ 2 + 0(\D E \ 3 ), (A.l) 

l CMB ^ ^ 

where is the cosmological kinematic signal as measured by 
a pencil beam, x is the pointing direction in the Ecliptic ref- 
erence frame and D E = (v Sun + Vpianck)/c in the same refer- 
ence frame. The equation accounts 13 for both the CMB dipole 
and the relativistic transverse Doppler effect, since the latter can 
be considered as a correction to the quadrupole proportional 
to 7cMB|vsun/c| 2 /2 ~ 19 //K. (We can safely discard -^\D E \ 2 , 
since it is a constant and plays no role in the calibration.) We 
denote with x\, = (x\,,y\,,Zb) the observing direction in the ref- 
erence frame of the beam. Then the convolved dipole signal is 



A7XQ 
To 



MB J 4k 



j 

J 4) 



N dQ.B(x h )x h -\J(t)D E 



~ f dnB(x h )(x b -U(t)D E f, 

1 J4x 



(A.2) 



13 Higher-order terms result in corrections at least a thousand times 
smaller and are therefore neglected. 



where U(f) is the matrix which transforms vectors from the 
Ecliptic reference frame to the beam reference frame, and N is a 
normalization constant: 



AT 



= f dQB(x h ). 

J4n 



(A3) 



The directional moments of the beam are defined as 

M/, m ,, 



„ = N f d£lB(x h )x' h yZz" h . (A.4) 

J4n 



They depend only on the shape of the beam and their value does 
not change during the mission. If we set 



S = (Mi i0 ,o,Mo,i,o,Mo,o,t) 
(a vector in the beam reference frame) and 



fM 2 n, , M, 



M 



1,0,1 



M u ,o, M ,2,o, M ,i,i 
Mi,o,i, M ,i,i, M ,o,2J 



(A.5) 



(A.6) 



(a symmetric matrix), then Eq. (A. 1) reduces to the bilinear form 
AT(t) 



'- CMB 



I* U:.<: S + i»;.U:/: AU .':/> . 



(A.7) 



A.2. Numerical calculation of the directional moments 



The computation of the components of 5 (Eq. A.5) and A 
(Eq. A.6) was based on the integration of separate GRASP maps 
for both sidelobes and main beam, plus intermediated beams 
computed for each radiometer at its nominal central frequency, 
i.e., 30 GHz, 44 GHz, and 70 GHz. In this context, sidelobes are 
defined as the part of the beam pattern outside 5° from the beam 
optical axis. The inner part of the beam pattern defines the inter- 
mediate beam and main beam maps. 

Each map takes into account the small differences in the 
optical path for each radiometer, given its different location in 
the focal-plane and polarization. For sidelobes, we used GRASP 
maps with a resolution of 0.5° along the azimuthal angle and 
0.5° along the zenith angle. For intermediate and main beams, 
the resolution is 0.5° (azimuth) and V (zenith). We integrated 
Eq. (A.4) separately for the intermediate beam and the sidelobes, 
and then we summed them together. (The same was done for the 
normalization constant N, Eq. A.3.) To compute the integral, we 
converted GRASP maps into HEALPix maps (Gorski et al. 2005) 
by bilinear interpolation and applied a simple midpoint algo- 
rithm. We chose the resolution of the HEALPix maps in order 
to have a resolution much smaller than that of the correspond- 
ing GRASP maps. Sidelobes were gridded onto an HEALPix map 
with A^ s ide = 1024, and intermediate beams onto a map with res- 
olution equivalent to A^ s id e = 65 536. 

In order to check the convergence, we computed the value of 
M/ jm> „ a number of times for increasing values of AT s ide- All the 
integrals showed a regular convergence curve. The largest A^de 
was fixed by having M, >mj „(A^ s ide)/M ; , m> „(A r S ide/2) - 1 < 10~ 3 for 
the moment with the slower convergence (including AO- 

The beam pattern changes over the bandpass of each ra- 
diometer, causing the moments M/ jmj „ to change slightly with fre- 
quency. However, the change mainly affects the sidelobes, which 
contribute to less than 0.5 % of the full-sky average. For temper- 
ature maps it is therefore enough to use a monochromatic cor- 
rection. 
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